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PREFACE 


The  various  Naval  laboratories,  under  sponsorship  of  the  Special 
Projects  Office,  have  undertaken  a  program  of  exploratory  scientific 
investigation  calculated  to  influence  the  future  Navy’s  at-sea  deterrence 
capability.  A  program  of  directed  research  of  some  duration  is  antici¬ 
pated  (and  required)  to  put  together  elements  which  are  fundamental  to 
the  formulation  of  a  future  sea-based  deterrence  weapon  system 
concept. 

An  appropriate  administrative  iramework  has  been  established  at 
NRL  to  sponsor,  on  a  Laboratory-wide  basis,  the  conduct  of  applicable 
e.xploratory  research.  Investigations,  matched  to  existing  NRL  scien¬ 
tific  programs  and  skills,  have  been  initiated  within  this  framework. 
They  are  broadly  directed  toward  a  consolidation  of  scientific  knowl¬ 
edge  and  technological  c  a  p  a  i  1  i  1 1  e  s  applicable  to  the  initiation  of  a 
firmly  scheduled  development  by  the  Navy  of  a  sea-based  deterrence 
weapon  systezn  matching  the  operational  requirements  of  its  day. 

This  report  is  one  of  a  series  intended  to  display  the  NRL 
research  activity  directed  to  support  this  follow-on,  sea-based  deter¬ 
rence  weapon  system  concept.  The  five  papers  included  in  this  report 
are:  A  Graphical-Numerical  Method  for  the  Transient  Response  of 
Nonlinear  Systems:  Influence  ot  Crystal  Surface  on  the  Optical  Trans¬ 
mission  of  LiF  in  the  Vacuum  Ultraviolet;  Small  Fatigue  Cracks  as 
Fracture  Origins  in  Tests  ol  High-Strength  Steel  Sheet:  Application  of 
Fracture  Mechanics  to  Adhesive  Joints:  and  A  New  Submarine  Detection 
System. 

Peter  Waterman 
Staff  Adviser  for  FBM 
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A  GRAPHICAL-NUMERICAL  METHOD  FOR  THE 
TRANSIENT  RESPONSE  OF  NONLINEAR  SYSTEMS 
[  Unclassified] 


P.  F.  Cunniff 

The  Naval  Research  Laboratory  (Mechanics  Divisionjhas 
been  conducting  transient-response  studies  of  nonlinear  one- 
and  two-degree-o£-freedom  nonautonomous  systems  by  a 
graphical-numerical  method,  which  is  a  modification  of  the 
general  phase-plane  method,  in  conjimction  with  numerical 
computations.  Present  analytical  procedures  are  adequate 
for  response  studies  of  linear  structures,  but  there  is  need 
for  new  and  better  procedures  for  analyzing  the  effects  of 
structural  nonlinearities.  This  is  a  step  in  the  continuing 
problem  of  the  evaluation  of  the  effects  on  equipment,  on 
ships  or  missiles,  when  subjected  to  the  various  shock  and 
vibration  environments  where  the  often-nonlinear  struc¬ 
tural  “feedback"  effects  must  be  recognized. 

. FBM  Staff 


INTRODUCTION 

One  of  the  more  important  decisions  which  the  engineer  must  make  in  the  field  of 
structural  dynamics  is  the  choice  of  model  to  serve  as  a  representation  of  a  structure  sub¬ 
jected  to  dynamic  forces.  One  model  configuration,  which  is  used  herein,  consists  of  con¬ 
centrated  masses  connected  in  a  straight  line  by  weightless  springs  and  dashpots.  The 
behavior  of  the  model  under  prescribed  forcing  functions,  which  are  called  inputs,  provides 
the  analyst  with  the  dynamic  characteristics  of  the  structure;  such  behavior  is  classified 
as  transient  response,  steady-state  vibrations,  and  random  vibrations.  Each  classification 
of  response  depends  on  the  input  to  the  system.  In  what  follows,  only  transient  response  is 
considered. 

The  physical  characteristics  of  the  models  as  represented  in  the  differential  equations 
of  motion  can  lead  to  either  linear  or  nonlinear  model  systems.  The  dynamic  analysis  of 
linear  elastic  structures  has  been  previously  reported  (1,2).  For  those  models  which  have 
nonlinear  characteristics  of  either  mass,  damping,  or  restoring  forces,  the  transient  solu¬ 
tion  of  the  differential  equations  of  motion  is  usually  unknown  as  a  combination  of  elemen¬ 
tary  functions,  and  approximate  techniques  are  resorted  to.  These  methods  are  generally 
numerical,  graphical,  and  electrical  analog  solutions.  Graphical  and  numerical  methods 
have  been  developed  for  single -degree -of-freedom  systems  and,  in  many  instances,  are 
applied  only  to  autonomous  systems. 

It  is  the  purpose  of  this  report*  to  present  a  reasonably  rapid,  yet  fairly  precise, 
technique  for  the  approximate  transient  solution  of  nonlinear  multi-degree-of-freedom 
systems.  The  nonlinear  components  of  the  systems  under  study  are  assumed  to  be  cubic 
hardening  and  cubic  softening  springs. 


*A  slightly  more  expanded  form  of  the  material  contained  in  this  report  is  available  as 
NRL  Report  6785,  of  the  same  title. 
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GRAPHICAL-NUMERICAL  METHOD 
General  Equations  of  Motion 

The  equations  of  motion  for  the  model  configuration  under  study  which  consists  of  r 
concentrated  masses  can  be  reduced  to 

Xj+it'?(X|  +  5>j)  =  0.  i-1.2 .  r  (1) 

where 

~  ’‘i  •  ’‘i'l’  *i+l'  ’‘i  ’  ’‘i-l- 

and  the  dots  refer  to  differentiation  with  respect  to  time.  A  graphical  solution  for  this 
system  can  be  obtained  on  r  phase-planes  by  extending  the  phase-plane  method  (3)  in  con¬ 
junction  with  simple  numerical  computations. 


Single  Degree  of  Freedom 

The  phase-plane  method  is  normally  applied  to  nonlinear  single-degree-of-freedom 
systems  and  has  been  extended  to  linear  two-degree -of -freedom  systems  by  superimposing 
normal  coordinates  (4).  The  differential  equation  of  motion  for  a  single-degree-of-freedom 
system  might  be  written  as 


X  +  (X  +  fj)  =  0.  (2) 

A  geometrical  interpretation  of  Eq.  (2)  can  be  obtained  by  introducing  new  coordinates 
which  are  called  the  phase-plane  coordinates.  Let  x  =  x  and  u  =  x/w.  Equation  (2)  then 
reduces  to  the  form 


The  5  has  the  same  dimension  as  x  and  may  be  represented  along  the  x-axis  in  the  phase- 
plane  as  shown  in  Fig.  1(a),  The  slope  of  the  normal  to  the  line  PQ  is  -u/fx  +  S).  Since 
this  slope  is  the  same  as  Eq.  (3),  the  phase  trajectory  at  P  will  be  normal  to  the  direction 
PQ. 


Fig.  1  -  The  phase-plane:  (a)  phase-plane 
coordinates  representing  the  differential 
equation  of  motion,  *  +  +  5)  =0,  where 

X  -  X  and  u  =  (b)  geometrical  represen- 

lation  of  the  phase  trajectory  dx/du.  From 
(a)  we  see  that  the  slope  of  the  normal  to 
PO  is  -u  /(  X  +  *  ),  This  is  also  given  by  Eq.  (3) 
asdx/du.  Thus  the  phase  trajectory  dx /du 
at  P  is  normal  to  the  direction  PQ. 
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Assume  that  s  is  held  at  its  mean  value  for  a  finite  time  interval  At.  Equation  (3) 
may  now  be  integrated  and  more  information  about  the  phase  trajectory  obtained.  Re¬ 
writing  Eq.  (3)  as 


(X  +  8)  dx  +  u  du  =  0 


and  integrating  leads  to 

(x  +  8)^  +  =  constant  =  (PQ)^,  (4) 

as  can  be  seen  in  Fig.  1(a).  According  to  Eq.  (4),  the  instantaneous  phase  trajectory  is  a 
circular  arc  passing  through  p  with  its  center  at  x  =  -8,  u  =  0,  i.e.,  at  Q. 

To  derive  a  relationship  between  the  circular  arc  segment,  beginning  at  P,  and  the 
time  increment  At,  note  that 


and 


dt 

ecu 


(5) 


Consider  Fig.  1(b),  which  shows  the  differential  geometry  of  the  phase  trajectory.  Let  the 
average  radius  of  curvature  of  the  step  be  given  as  qp  =  p.  Then, 

dO  =  —  =  /dx^  +  du^ 


This  leads  to 


^  l/u^  +  (x  +  8)2 
./i 

-  ^  /  ^dx^  _  ^ 


by  means  of  Eq.  (3).  Substitution  of  this  last  equation  into  Eq.  (5)  yields 


dt  =  ‘5®. 


Since  the  mean  value  of  8  was  used  in  the  derivation,  this  may  be  written  as 

AS  -  a;  At.  ^  (6) 

Positive  time  is  therefore  represented  by  a  counterclockwise  angular  variation  of  the  nor¬ 
mal  PQ.  If  a  linear  or  nonlinear  second-order  ordinary  differential  equation  is  reduced  to 
this  delta  form,  it  can  be  integrated  graphically  by  construction  of  successive  circular  arc 
segments  whose  centers  are  located  on  the  x-axis  at  various  step  values  of  x  =  -8.  The 
locus  of  these  circular  arc  segments  is  the  approximate  phase  trajectory  and,  therefore, 
the  desired  solution  of  the  differential  equation. 


Method  of  Construction 

The  construction  procedure  of  the  phase  trajectory  by  the  delta  phase -plane  method 
for  nonautonomous  single-degree-of-freedom  systems  is  now  described.  As  an  example. 
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consider  a  simple  damped  oscillator  with  a  cubic  hardening  spring  and  a  prescribed  foun¬ 
dation  velocity  input  (y„  ).  The  equation  of  motion  is 

mic  +  cx  +  kx  +  /3x^  - 


This  is  written  in  the  S  form  as 


X  +  du2(x  +  8)  -  0. 

where 

8  =  -y2x3  +  2au  +  Zi  . 


The  constants  y  and  a  are  given  by  y  =  and  a  =  Assume  that  the  graphical 

solution  has  been  extended  to  step  n,  i.e.,  x„  and  u„  are  known  and  x„<,i  and  are 
desired.  Using  the  average  values  of  the  phase  coordinates  during  the  next  increment  of 
time,  the  8,,  equation  is  ‘ 


8 


n 


‘"(“n 


(7) 


Since  the  input  is  known  as  a  function  of  the  foundation  velocity,  the  foundation  accelera¬ 
tion  may  be  approximated  over  the  finite  time  interval  At  by  the  theorem  of  the  mean  of 
differential  calculus.  Thus, 


(Vo), 


At 


(8) 


So,  all  terms  in  the  right-hand  side  of  Eq.  (7)  are  known  except  and  u„+i  . 

As  a  first  trial  for  finding  and  u„n,  the  known  8n_i  from  the  previous  step  of 
the  solution  might  be  used  as  the  first  center  of  arc.  By  laying  off  the  prescribed  circu¬ 
lar  arc  A”  =  f.At  from  the  coordinate  (ii„,  x„)  by  means  of  a  compass  and  protractor,  pre¬ 
liminary  values  of  and  u  are  found  quickly.  Calculate  the  first  trial  value  of  s,, 
from  Eq.  (7).  Draw  a  new  circular  arc  from  (u^,  with  the  new  center  at  x  =  -8^  so  that 
a  new  coordinate  x^^j)  is  found.  Again,  use  this  coordinate  to  calculate  a  new 

This  iteration  process  continues  until  the  values  of  8„  converge.  This  completes  one  step 
of  the  solution. 

The  initial  conditions  for  a  system  are  known,  but  there  is  no  preceding  value  of  8 
established  at  the  beginning  of  a  problem.  To  calculate  a  trial  value  of  the  magnitude 
of  the  phase  coordinates  (uj.xj)  at  the  end  of  the  first  step  might  be  assumed  or  the 
initial  values  of  the  system  might  be  used  as  the  average  values  over  the  interval,  or  a 
Maclaurin  expansion  might  be  used  to  calculate  the  phase  coordinate  (uj,  xj)from  the 
initial  values. 


Two  Degrees  of  Freedom 

The  equations  of  motion  for  a  tv'o-degree-of-freedom  system  are  written  in  the  8 
form  as  follows: 


Xj+<aJ(Xi+Si)=0 

Xj  +  a>^(X2  +  82)  ~  0, 


and 


(9) 

(10) 
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where 

and 


Sj  Xj,  Xj, 

Xj,  t  ) 

Sjfxj.  Xj.  xj. 

*2-  t  ) 

For  example,  if  the  equations  of  motion  of  a  linear  undamped  system  are  reduced  in  terms 
of  the  relative  displacements,  the  s  expressions  are 


and 


Si 


+ 


Yo 


CO 


2 

1 


X 


1- 


The  graphical  integration  of  the  coupled  Eqs.  (9)  and  (10)  is  obtained  by  extending  the 
delta  phase-plane  method  to  two  phase-planes,  (u^,  xj)  and  (u^,  x^).  The  computational 
method  described  in  the  preceding  section  is  also  used.  If  the  solution  has  been  extended 
to  step  n,  (5i)„and  (  62)n  are  found  for  step  (n+i)  by  using  the  known  (Si)„_iand  (5  2)n-i^s 
the  initial  centers  of  circular  arcs.  This  gives  trial  values  of  [(xj)^^.„  (ui)„+,]and  [(xjln+j, 
and  (uj )„+,],  which  allows  (s^)^  and  (s^)^  to  be  computed.  The  procedure  is  repeated,  as" 
in  the  case  of  the  single-degree-of -freedom  system,  until  (&!)„  and  (^2)^  converge  to  their 
respective  values. 


EQUIVALENT  NUMERICAL  METHOD 

The  differential  equation  of  motion  for  a  linear  undamped  oscillator  subjected  to  a 
foundation  motion  is 


X  +  .V  2  X  =  -  y^. 

A  numerical  integration  technique  has  been  developed  (5)  for  the  solution  of  this  equation 
where  the  input  is  known.  For  a  finite  increment  of  time  .'.t  =  ii,  assume  that  the  initial 
conditions  for  the  (n  +  ij  interval  are  and  u^.  The  solution  at  the  end  of  the  interval  is 
then 


and 


’‘n+l  "  ’‘p. 


r*'  r 

sin  -h  -  i  1  sin  a.(h  -  T)  dT 


"  Xj,  sin  j.h  +  Uj,  cos 


“"-JJ  [m 


(T)  cos  c.(h  -  T)  dT. 


K  the  foundation  acceleration  is  defined  by  Eq.  (8),  the  approximate  solution  equations 
become 

(yo)„ 

+  l  ""  ’‘n  +  u„  Sin  coh  -  - —  (  1  -  CO  s  coh) 


■'n+1 


"  x„  sin  c*jh  +  u„  cos  coh  - 


(yo)„ 


r.2 


sin  <i.h. 


This  same  approach  is  now  used  for  the  differential  equation  of  motion  in  the  standard 
0  form,  that  is,  for 

X  +  co^X  =  -co2  h.  (11) 
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Assuming  that  5  is  held  constant  during  a  finite  time  increment,  the  solution  equations 
for  X  and  u  at  the  (n+  1)  interval  are 


and 


=  Xj^  cos  ii)h  +  Uj,  sin  a>h  -  *^^(1  -  cos  ojh) 
~  Xn  sin  ojh  +  COS  coh  “  5^  sin  a>h. 


(12) 

(13) 


The  first  two  terms  on  the  right-hand  side  of  Eqs.  (12)  and  (13)  are  known,  as  well 
as  the  constant  coefficients  of  the  third  term  The  procedure  consists  in  using 
from  the  previous  step  for  in  Eqs.  (12)  and  (13).  These  first  trial  values  of  Xn+^and 
are  used  to  find  a  trial  value  of  Sn.  Equations  (12)  and  (13)  then  give  new  values  of 
Xj^^j  and  which  are  used  to  compute  a  new  value  of  The  numerical  procedure  con¬ 
tinues  until  the  values  of  converge. 

To  show  that  this  numerical  procedure  is  the  equivalent  procedure  of  the  S  graphical 
construction  method,  rearrange  Eqs.  (12)  and  (13)  in  the  following  form: 


and 


Xj,4.j  "  x^  cos  ii)h  +  sin  cos  cijh 

<i)h'+  cos  cuh  -  sin  oih. 


(14) 

(15) 


Square  each  side  of  Eqs.  (14)  and  (15)  and  add  the  two  resulting  equations.  Upon  simpli 
fying,  the  equation  reduces  to 


«n)' 


=  (’‘n  +  *n)' 


(16) 


This  is  precisely  the  equation  of  a  circle  whose  center  is  at  =  -  5-,  u^+j=  0,  and  whose 
radius  squared  is  (x„  +  s^)^  +  u^.  It  is  concluded  that  Eqs.  (12)  and  (13)  give  the  same 
result  as  the  previously  described  graphical -numerical  method  without  a  drawing  error. 

If  Eqs.  (9)  and  (10)  are  rearranged  with  the  5  terms  on  the  right-hand  side  of  the 
equations,  the  solution  at  the  (n+l)  interval  is  established  in  a  similar  manner  as  in  the 
case  of  a  single -degree -of -freedom  system.  Thus, 


and 


(Xj)  -  (Xj)  cos  <iijh  +  (Uj)  sin  wjh  -  (1  "  cos  c.,jh), 

(Ui)n+i  ^  wjh  +  {“i)„  cos  wjh  -  (Sj)^  sin  oijh, 

(Xj)  =  C"2)  cos  c.;2h  +  (u  2 )  sin  cijjh  ~  (82)^  (1  “  cos  cujh). 


(U2)  - (*2)  sin  oijh  +  ^2^  ~  (^2\  ‘^2^’ 

Once  again,  all  terms  in  the  right-hand  side  of  these  four  equations  are  known,  except 
(a,)  and  (a,)  where 

An  ^  n 


and 


"  El  [(Xl)„+i. 

(C2 

{Uj 

1 - 1 

4-1 

C 

The  same  iteration  process  is  followed  for  the  convergence  to  (a^)  and  (a 2)  as  in  the 
case  of  in  the  single -degree -of -freedom  system. 
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RESULTS 

For  a  mechanical  model  which  idealizes  a  nonlinear  system  under  a  dynamic  input, 
it  is  generally  important  to  know  the  relative  displacements  during  transient  motion. 
Example  problems  are  presented  for  one-  and  two -degree -of -freedom  systems  which 
apply  the  methods  developed  in  this  report  for  finding  the  relative  displacements. 

The  graphical -numerical  method  is  applied  to  the  four  systems  shown  in  Fig.  2.  In 
each  case  the  system  is  initially  at  rest  and  the  input  is  the  foundation  velocity  shown  in 
Fig.  3.  The  time  duration  of  the  input  is  t  =  0.500  second,  and  the  solution  curves  for  each 
example  extend  to  t  =  1.000  second. 


Fig.  2  -  Four  possible  systems  to  which  the  phase-plane  method 
is  applied  in  this  report:  (a)  linear,  single-dcgree-of-freedom 
system  with  viscous  damping;  (b)  nonlinear,  single-degree-ol- 
freedom  system  with  both  viscous  damping  and  a  cubic  hardening 
spring;  (c)  nonlinear,  two-degree-of-freedom  system  without 
damping  and  with  a  cubic  hardening  spring  between  mj  and  m,  ;  and 
(d)  nonlinear,  two-degree-of-freedom  system  with  viscous  damp¬ 
ing  and  a  cubic  hardening  spring  between  m,  and  mj  . 


Fig.  3  -  Example  of  a  prescribed  foundation 
velocity  input  (The  time  duration  of  the  input 
for  this  example  is  0.500  sec.)  For  a  simple 
damped  oscillator  with  a  cubic  hardening  spring 
and  a  prescribed  foundation  velocity  input,  the 
equation  of  motion  is  mx+  cx  +  fcx  +  4x^  =  -  my^. 


8 


NAVAL  RESEARCH  LABORATORY 


The  construction  of  the  phase -plane  trajectories  by  the  graphical -numerical  method 
is  obtained  by  using  a  compass  to  lay  off  the  required  circular  arc  and  a  protractor  to 
measure  the  angle  of  increment  A0.  An  engineer’s  scale  measures  the  coordinate  points 
on  the  phase-plane.  The  solutions  for  the  nonlinear  systems  are  compared  with  their 
respective  approximate  solutions  on  an  analog  computer  at  NRL.  The  computer  results 
for  each  problem  arc  in  the  form  of  plots  of  relative  displacement  x  versus  time  t.  The 
graphical -numerical  solution  and  the  exact  solution  of  the  linear  system  shown  in  Fig.  2(a) 
are  presented  in  the  Appendix  and  are  compared  with  the  approximate  results  of  the  analog 
computer  for  this  problem  (Fig.  Al(b),  Appendix). 

The  system  shown  in  Fig.  2(b)  is  solved  in  Example  1(b)  by  the  numerical  equivalent 
procedure  of  the  graphical -numerical  method.  This  provides  a  check  on  the  drawing  error 
for  this  problem. 


Example  1(a) 

Fig.  2(b)  is  a  model  of  a  nonlinear  single-degree-of-freedom  system  with  viscous 
damping  and  a  cubic  hardening  spring.  The  S  equation  is 

5  =  2au  +2L 

2  2 

The  parameter  values  are 

-  n  Ih  -  scc^  ,  _  »  lb . 

m  -  0.20  - k  -  77.2  — : 

in.  in. 


So, 


0.7S  ,3  =  15.4 


b  =  0.191  n  +  0.199 

.  2 


An  increment  of  time  At  =  0.025  second  was  used  throughout  the  solution  so  that  A6'=  28.13 
degrees.  The  equation  is 


2  8  9.650 

Since  the  mass  is  initially  at  rest,  the  initial-time  coordinate  conditions  are 


=  0  and  ii^,  -  “y^  w  -  -1.272  in. 

Figure  4(a)  is  the  phase-plane  trajectory  composed  of  40  steps  in  the  graphical- 
numerical  solution.  As  a  representative  step  in  the  solution,  consider  the  construction 
of  the  arc  segment  from  point  6  to  point  7.  At  the  end  of  step  6,  =  0.391  inch,  Ug  =  0.826 

inch,  and  S5  =  -  0.006  inch.  Use  as  a  first  trial  to  find  Xy  and  U7.  This  construction  is 
shown  in  Fig.  4(b),  and  xy  ~  0.737  inch  and  u,  -  0.545  inch.  Using  these  values,  Sg  is 
calculated  to  be  0.026  inch.  A  new  arc  is  drawn  as  shown  in  Fig.  4(b)  so  that  xy  -  0.731 
inch  and  uy  -  0.532  inch.  The  second  value  of  Sg  is  0.026  inch.  Since  the  engineer’s  scale 
can  be  read  to  within  approximately  ±0.002  inch,  the  construction  of  point  7  is  complete. 

Figure  4(c)  is  the  analog  computer  response  curve.  As  a  comparison  between  the  two 
solutions,  the  points  plotted  in  Fig.  4(c)  represent  the  x  coordinates  of  Fig.  4(a)  at  0.025- 
second  intervEds.  Solution  points  obtained  by  the  numerical  equivalent  method  as  presented 
in  Example  1(b)  are  shown  in  Fig.  4(c). 


X  (IN.) 
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>  (IN.) 


»  (IN.) 


Fig.  4  -  Solution  of  the  nonlinear,  single-degree-of-freedom  system  with 
viscous  damping  and  a  cubic  hardening  spring,  (a)  phase-plane  trajectory 
consisting  of  40  steps  in  the  graphical-numerical  solution  with  As  =  28.13 
degrees;  (b)  illustration  of  two  trials  in  the  evaluation  of  the  step  from  (x^, 
Uj)  to  (x,,u,);  (c)  analog-computer-approximation  response  curve  show¬ 
ing  values  of  x  as  functions  of  t.  Selected  results  from  the  graphical- 
numerical  method  are  shown  for  0.025-second  intervals  as  circles,  the 
triangles  representing  some  selected  exact  numerical  results. 
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Example  1(b) 


The  numerical  equivalent  procedure  of  the  graphical -numerical  method  is  used  to 
solve  Example  1(a).  Since  At  -  h  =  0.025  second  throughout  the  problem,  the  solution 
equations  are 


and 

where 


’‘n+l  “  0.88178  x„  +  0.47166  "  0.11822 


""  -0.47166  x„  +  0.88178  -  0.47166  5^, 


=  0.09544  +  u„)  +  0.02494 


s  a. 

^  +  = - 

n '  /\ 


9.650 


The  results  of  this  method  are  plotted  in  Fig.  4(c)  only  where  there  is  a  noticeable  differ¬ 
ence  from  the  values  of  x  in  Example  1(a). 


Example  2 


Figure  2(c)  is  a  model  of  a  nonlinear,  two -degree -of -freedom  system  without  damp¬ 
ing  and  with  a  cubic  hardening  spring  between  masses  mj  and  The  i  equations  are 

ko 


"2  ..3  .  yo 


and 


^  k.  2  kj  2  2 

o-i 


I  _  "'2  .  /  1  .  1  \  ^’2  3 

2  ftij  2  ^n-j  2  ^  2  ' 


The  parameter  values  are 


=  0.20  n,.  =  0.10  ibjLHs! 

*  I  r\  ^  I  r» 


Hence, 


and 


k  =  77,2  kj  -  38.6  =  7.7 

in.  in.  ^  jp,  3 

V 

5j  =  -  0.500  Xj  -  0-100  ^2  *  — 

^2  =  0.500  Xj  +  0.299  Xj  -  Xj. 


Since  =  0,2  -  19.65  rad/sec,  the  angle  of  increment  is  the  same  on  each  phase-plane, 
i.e.,  A^i  =  A^j  =  A(5.  Figures  5(a)  and  5(b)  are  phase  trajectories  for  A0  =  14.07  degrees 
which  required  80  steps  in  the  solution.  Figures  5(c)  and  5(d)  are  the  response  curves  of 
the  analog  computer.  The  values  of  (xi)„  and  (x2)„  from  Figs.  5(a)  and  5(b)  are  plotted  on 
Figs.  5(c)  and  5(d)  at  time  intervals  At  =  0.025  second  as  a  comparison  of  the  solutions. 
Response  values  for  A6'  =  28.13  degrees  by  the  graphical -numerical  method  are  also  shown 
in  the  figures. 


Example  3 

Figure  2(d)  is  a  model  of  a  nonlinear  two -degree -of -freedom  system  with  viscous 
damping.  The  problem  is  the  same  as  Example  2  with  the  addition  of  the  damping  constants, 
Cj  =  0.75  Ib-sec/in.  and  C2  =  0.25  Ib-sec/in.  The  8  equations  are 

y 

Sj  =  "0.500  Xj  “  0.100  j  +  0.191  Uj  -  0.064 


(IN  ) 
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Fig.  5  -  Solution  of  the  nonlinear,  two-degree-of-£reedom  system  without 
damping  and  with  a  cubic  hardening  spring,  (a)  phase-plane  trajectory 
consisting  of  80  steps  in  the  graphical  numerical  solution  for  the  evalua¬ 
tion  of  jt  and  ij  for  mj  with  —  14,0V  degrees^  (b)  the  same  as  (a),  except 
the  values  are  for  nij;  (c)  analog-computer-approximation  response  curve 
showing  values  of  x  as  functions  of  t ;  (d)  the  same  as  (c),  except  the  curve 
is  for  values  of  x  for  the  mass  mj.  In  (c)  and  (d)  the  circles  are  selected 
results  obtained  by  the  graphical-numerical  method  at  0.02  5- second  inter¬ 
vals  for  he  -  28.13  degrees,  and  the  triangles  are  some  selected  exact 
numerical  results  obtained  for  he  =  14,07  degrees. 
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(d) 

Fig.  5  (cont’d)  -  Solution  of  the  nonlinear,  two-degree-of-freedom  system 
without  damping  and  with  a  cubic  hardening  spring,  (a)  phase-plane  trajec¬ 
tory  consisting  of  80  steps  in  the  graphical  numerical  solution  for  the  evalua¬ 
tion  of  *  and  u  for  m,  with  =  14.07  degrees;  (b)  the  same  as  (a),  except 
the  values  are  for  mj;  (c)  analog-computer-approximation  response  curve 
showing  values  of  x  as  functions  of  t ;  (d)  the  same  as  (c),  except  the  curve 
is  for  values  of  x  for  the  mass  mj.  In  (c)  and  (d)  the  circles  are  selected 
results  obtained  by  the  graphical-numerical  method  at  0.025- second  inter¬ 
vals  for  Afl  =  28.13  degrees,  and  the  triangles  are  some  selected  exact 
numerical  results  obtained  for  AS  =  14,07  degrees. 


and 

^2  "  0.500  Xj  0-299  *1  “  *1  +  0-191  Uj  “  0.191  u^. 

Figures  6(a)  and  6(b)  are  the  graphical  solutions  consisting  of  80  steps  at  M  =  14.07 
degrees,  while  Figs.  6(c)  and  6(d),  are  the  solution  curves  obtained  from  the  analog  com¬ 
puter.  The  values  of  (xi)„  and  (x2)„  from  Figs.  6(a)  and  6(b)  are  plotted  at  time  intervals 
At  =  0.025  second  on  Figs".  6(c)  and  6(d)  for  comparison  of  the  solutions.  Response  values 
for  A6  =  28.13  degrees  by  the  graphical -numerical  method  are  also  shown  in  the  figures. 


DISCUSSION 

The  construction  of  the  phase-plane  trajectories  by  the  graphical -numerical  method 
required  a  protractor,  an  engineer’c  scale,  a  set  of  triangles,  and  a  compass,  while  a 
desk  calculator  was  used  for  computing  the  numerical  values  of  the  S’s.  A  scale  of 
1  inch  =  0.400  inch  was  selected  for  each  phase-plane  coordinate  system  so  that  the  inves¬ 
tigator  was  able  to  read  a  coordinate  to  within  0.002  inch.  (The  construction  of  a  circular 
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Fig.  6  -  Solution  of  the  nonlinear,  two-degree-of-freedom  system  with 
viscous  damping  and  a  cubic  hardening  spring,  (a)  phase- plane  trajectory 
consisting  of  80  steps  in  the  graphical -numerical  solution  for  the  evaluation 
of  X  and  u  for  mj  with  Ae  =  14.07  degrees;  (h)  the  same  as  (a),  except  the 
values  are  for  (c)  analog-computer-approximation  response  curve 
showing  values  of  x  as  functions  of  t;  (d)  the  same  as  (c),  except  the  curve 
is  for  values  of  x  for  mass  mj.  In  (c)  and  (d)  the  circles  are  selected 
results  obtained  by  the  graphical-numerical  method  at  0.025-second  inter¬ 
vals  for  Ae  =  28.13  degrees,  and  the  triangles  are  some  selected  exact 
numerical  results  obtained  for  Ae  =  14.07  degrees. 
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Fig.  6  (cont’d)  -  Solution  of  the  nonlinear,  two-degree-of-freedom  system  with 
viscous  damping  and  a  cubic  hardening  spring,  (a)  phase-plane  trajectory 
consisting  of  80  steps  in  the  graphical-numerical  solution  for  the  evaluation 
of  X  and  u  for  r- ,  with  A3  =  14.07  degrees;  (b)  the  same  as  (a),  except  the 
values  are  for  m,;  (c )  analog-computer-approximation  response  curve 
showing  values  of  x  as  functions  of  t  ;  (d)  the  same  as  (c),  except  the  curve 
is  for  values  of  x  for  mass  inj.  In  (c)  and  (d)  the  circles  are  selected 
results  obtained  by  the  graphical-numerical  method  at  0.025- second  inter¬ 
vals  for  A3  =  28.13  degrees,  and  the  triangles  are  some  selected  exact 
numerical  results  obtained  for  A3  =  14,07  degrees. 


arc  segment  for  each  step  required  that  the  trial  values  of  the  3  ’s  converge  to  within 
±0.002  inch.)  The  protractor  which  measured  the  angle  of  increment  is  calibrated  to 
0.5  degree.  It  appears  that  an  angle  can  be  read  approximately  to  within  0.2  degree. 

The  example  of  a  linear  single-degree-of -freedom  system  in  the  Appendix  serves 
two  purposes.  The  exact  values  of  x  are  compared  with  the  results  of  the  graphical- 
numerical  method  (Fig.  Al(a))  and  the  analog  computer  (Fig.  Al(b)).  This  provides  a  check 
on  the  accuracy  of  the  analog  computer.  The  values  of  relative  displacement  using  the 
graphical-numerical  method  are  plotted  at  increments  of  0.025  second  on  Fig.  Al(b).  The 
exact  values  of  x  are  computed  at  intervals  of  0.050  second  and  are  plotted  on  the  figure 
where  there  is  a  noticeable  difference  from  the  graphical -numerical  results.  The  greatest 
differences  between  the  exact  solution  and  the  computer  curve  occur  at  t  =  0.500  second 
and  t  =  0.750  second;  this  difference  is  approximately  0.033  inch.  Otherwise,  the  values 
compare  quite  favorably. 

The  results  of  the  nonlinear  systems  of  Examples  1,  2,  and  3  are  shown  in  Figs.  4,  5, 
and  6.  In  each  case  the  graphical -numerical  results  agree  favorably  with  the  analog  com¬ 
puter  curves.  The  iteration  process  for  each  example  generally  required  two  trials  in 
each  step  for  the  5’s  to  converge  to  within  ±0.002  inch. 

The  numerical  equivalent  procedure  was  used  in  Example  1(b)  to  find  the  drawing  error 
in  the  graphical -numerical  method  applied  to  the  system  of  Example  1(a).  The  results  show 
that  the  two  methods  are  in  agreement,  as  seen  in  Fig.  4(c). 
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CONCLUSIONS 

If  the  equations  of  motion  for  a  nonlinear  system  are  reduced  to  the  standard  delta 
form,  as  in  Examples  1(a),  2,  and  3,  the  graphical-numerical  method  is  readily  adaptable 
for  finding  the  approximate  displacement  for  a  given  transient  input.  For  nonautonomous 
systems  whose  input  is  in  the  form  of  a  foundation  velocity,  the  corresponding  foundation 
acceleration  at  finite  increments  of  time  can  be  approximated  by  the  theorem  of  the  mean 
of  differential  calculus.  The  graphical -numerical  method  is  straightforward  to  apply, 
and  the  time  required  to  solve  a  problem  depends  upon  the  number  of  variable  terms  in 
the  5  equation  and  the  degree  of  accuracy  desired.  From  the  experience  gained  in  the 
example  problems,  the  proposed  technique  of  using  from  the  previous  step  as  the 
first  trial  of  has  been  useful  for  obtaining  a  rapid  convergence  of  The  following 
listing  shows  the  approximate  time  required  for  one  person  to  solve  each  of  the  first 
three  examples. 


Example 

— 

Angular  Increment 
(deg) 

Steps 

Time 

(hr) 

1(a) 

28.13 

40 

6-1/2 

2 

14.07 

80 

16 

3 

14.07 

80 

20 

The  numerical  equivalent  procedure  is  applied  with  the  aid  of  a  desk  calculator. 

This  procedure  has  shown  that  the  drawing  error  of  the  graphical-numerical  method  can 
be  quite  small,  as  indicated  in  Example  1(b).  For  problems  where  a  high  degree  of  accu¬ 
racy  is  required,  one  might  first  draw  the  phase  trajectories  for  one  or  two  large  constant 
angular  increments  as  a  preliminary  approach  to  the  solution.  With  this  information,  the 
numerical  procedure  for  a  small  angular  increment  could  be  applied  to  find  a  more  accu¬ 
rate  solution.  The  preliminary  graphical  trajectories  will  show  where  the  solution  curves 
have  sharp  changes  of  curvature  and  small  radii  of  curvature.  At  such  regions,  the  angu¬ 
lar  increment  of  the  numerical  solution  might  be  further  reduced  in  magnitude  for  greater 
accuracy. 
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APPENDIX 


SOLUTION  OF  LINEAR,  SINGLE-DEGREE -OF-FREEDOM  SYSTEM 


Figure  2(a)  is  a  model  of  a  linear,  single -degree -of -freedom  system  with  viscous 
damping.  The  s  equation  for  this  system  is 

S  =  2au  +  ~  . 

The  parameter  values  are 

mg  =  60  lbs;  k  =  60  Ib/in.  :  c  =  0.3  . 

in. 

A  constant  increment  of  time,  At  =  0,025  second,  was  used  so  that  the  angle  of  arc  is 
A 5  i  28.13  degrees  for  each  step  of  the  solution.  Since  the  mass  is  initially  at  rest,  the 
initial -time  coordinate  conditions  are 


Xg  =  0  and  =  -y^/o.  =  -  1.272  in. 

Figure  Al(a)  is  the  phase-plane  trajectory  composed  of  40  steps  in  the  graphical- 
numerical  solution,  while  Fig.  Al(b)  is  the  response  solution  from  the  analog  computer. 
The  mathematical  solution  for  this  linear  problem  in  terms  of  the  foundation  velocity  is 

X  =  x(o)  cos  pt  *  — ^ —  [aaix{o)  +  x(o)]  sin  pt 


sin  pt  + 


r 


y  /j)  e-<iw(t-T> 


—  sin  p{  t  -T)  “  cos  p(  t  -T) 


dT. 


Values  of  X  at  time  increments  of  0.050  second  were  calculated  with  this  equation  and  are 
plotted  on  Fig.  Al(b)  where  there  is  a  noticeable  difference  from  the  values  of  the  graphical- 
numerical  method. 
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Fig.  A1  -  Solution  of  the  linear,  single-degree-of-freedom  system  with 
viscous  damping,  (a)  phase-plane  trajectory  consisting  of  40  steps  in 
the  graphical-numerical  solution  with  iff  =  28.13  degrees;  (b)  analog- 
computer-approximation  response  curve  showing  values  of  x  as  func¬ 
tions  of  t.  The  circles  on  (b)  represent  selected  results  obtained  by 
the  graphical-numerical  method  for  =  28.13  degrees  and  for  0.025- 
second  intervals,  and  the  triangles  are  some  selected  exact  numerical 
results  for  0.050-second  intervals. 


SYMBOLS 


cj  -  viscous  damping  constant  between  m.  andm^., 
fjQ  -  noniinear  spring  force  between  m  and  the  foundation: 

k^Xj  +  PlXj 

fji  -  nonlinear  spring  force  between  mj  and  m^;  k2X2  +  P2X2 
h  -  finite  time  increment:  At 
k.  -  spring  constant  between  idj  and  m._j 
rtij  -  ith  mass 

p  ”  oj/ 1  - 

t  -  independent  time  variable 
u-  “  scaled  system  velocity;  Xj/«j 

-  relative  displacement  between  and  mj_i 
y .  -  absolute  displacement  of  m. 

-  absolute  displacement  of  the  foundation 

“i  ~  2tlljr<,.  ) 

-  coefficient  of  nonlinear  term  of  spring 


-  phase-plane  delta  for  the  ith  equation  of  motion 
A-.:  -  angular  increment  of  circular  arc  of  the  phase-plane 

a,?  -  kj  / m. 
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INFLUENCE  OF  CRYSTAL  SURFACE  ON  THE  OPTICAL 
TRANSMISSION  OF  LiF  IN  THE  VACUUM  ULTRAVIOLET 

[  Unclassified] 

D,  A.  Patterson  and  W.  H.  Vaughan 

The  Naval  Research  Laboratory  (Solid  State  Divisiun) 
has  been  concerned  with  the  study  of  imperfections  and  the 
effects  of  imperfections  on  the  properties  of  solids.  These 
studies  have  emphasized  optical,  etch,  and  x-ray  techniques 
for  gross  imperfections.  A  sensitive  light  scattering  appa¬ 
ratus  has  been  built  and  is  being  used  to  study  impurities  in 
the  alkali  halides.  A  second  very  useful  technique  is  a  sen¬ 
sitivehardness  test  which  has  been  used  to  study  the  motion 
of  dislocations.  The  materials  being  investigated  include  a 
variety  of  water-soluble  crystals  ranging  from  simple  cubic 
alkali  halides,  both  pure  and  containing  various  additives, 
and  complex  structures  such  as  ammonium  phosphate. 

. FBM  Staff 


INTRODUCTION 

The  optical  properties  of  lithium  fluoride  single  crystals  have  been  studied  for  more 
than  fifty  years,  most  recently  by  Kato  (1)  et  al.  Since  those  crystals  are  transparent 
farther  into  the  vacuum  ultraviolet  (to  about  1050A)  than  any  other  material,  they  have 
found  a  great  many  uses  in  that  spectral  region  as  windows  for  light  sources  and  detec¬ 
tors,  substrates  for  study  of  thin  films,  and  as  prisms  and  lenses.  Since  the  vacuum 
ultraviolet  transmission  of  LiF  is  substantially  reduced  by  mechanical  polishing,  thin 
optical  plates  are  usually  obtained  by  cleavage.  Hence,  the  mechanical  properties  of  the 
crystal  have  an  important  effect  on  the  ease  witJi  which  optical  plates  are  obtainable. 

Some  LiF  crystals  recently  received  from  the  Harshaw  Chemical  Company  have  been 
found  to  be  much  more  difficult  to  cleave,  and  the  present  invc-stigation  started  out  simply 
to  improve  the  cleavage  properties  by  addition  of  divalent  impurities,  cleaving  at  low  tem¬ 
perature,  exposure  to  ionizing  radiation,  or  by  some  combination  of  the  three.  However, 
inconsistencies  in  the  optical  transmission  measurements  on  successive  days  had  to  be 
accounted  for  and  became  the  principal  subject  of  this  paper. 


EXPERIMENTAL 

The  optical  measurements  in  the  vacuum  ultraviolet  were  made  using  a  half-meter 
grating  spectrometer  of  the  Seya-Namioka  type  manufactured  by  the  McPherson  Instru¬ 
ment  Corporation.  The  light  source  was  a  dc  discharge  with  hydrogen  flowing  continuously 
through  the  capillary  into  the  vacuum  chamber.  Optical  transmission  values  were  com¬ 
puted  at  the  peaks  of  strong  lines  to  minimize  corrections. 

Samples  from  a  number  of  different  LiF  crystals  were  studied,  but  those  with  which 
the  data  in  this  paper  were  obtained  are  all  from  the  same  block  unless  otherwise  stated. 
They  were  sliced  from  2  x  1/2  x  1/2  inch  crystal,  purchased  in  1958  from  Harshaw  Chem¬ 
ical  Company,  which  cleaved  very  easily  into  good,  thin  plates.  The  samples  varied  in 
thickness  from  0.20  to  0.125  cm,  with  little  difference  being  observed  either  in  initial 
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properties  or  in  subsequent  behavior.  They  were  mounted  with  dots  of  light  vacuum  grease 
on  six -position  slides  which  could  be  translated  into  the  light  beam  for  measurement. 

High-humidity  exposures  were  made  in  a  small  (48  cu  in.)  plastic  box  half  filled  with 
water  above  which  the  sample  slide  was  supported.  Low-humidity  exposures  were  made 
in  a  desiccator  containing  PjOj.  These  will  be  referred  to  as  the  “wet”  and  “dry”  boxes, 
respectively.  Vacuum  treatments  were  done  in  a  vacuum  desiccator  (without  desiccant) 
with  a  liquid  nitrogen  cold  trap  in  the  pump  line.  Pressures  in  the  vacuum  treatments 
were  10"®  mm  Hg,  or  less,  and  will  be  referred  to  as  “good”  vacuum. 


RESULTS 

The  vacuum  ultraviolet  transmission  of  all  the  LiF  crystals  studied  has  been  found  to 
decrease  after  cleavage.  In  Fig.  1,  curve  (a)  shows  the  transmission  spectrum  of  a  typi¬ 
cal  crystal  immediately  after  cleavage.  The  same  crystal  97  days  later,  shown  in  curve 
(b),  has  suffered  a  fairly  general  loss  in  transmission  below  160  nifi.  This  loss  does  not 
appear  to  be  associated  with  the  growth  of  any  specific  absorption  band,  such  as  that  due 
to  OH"  shown,  for  example,  in  curve  (c)  in  a  different  crystal. 


Fig.  1  -  Optical  transmission  of  LiF:  (a) 
freshly  cleaved  crystal;  (b)  the  same  crys¬ 
tal  after  97  days  in  air;  and  (c)  another 
crystal,  showing  OH"  absorption  band  at 
136  (1360A). 


Since  no  specific  absorption  bands  appeared  as  the  transmission  decreased,  the  strong 
hydrogen  Lyman-a  line  at  1216A  was  chosen  as  a  convenient  point  at  which  to  follow  the 
decrease  as  a  function  of  time  and  treatment.  Figure  2  shows  a  typical  curve  of  trans¬ 
mission  (at  1216A)  versus  time  in  air.  The  apparent  “saturation”  (in  this  case  beyond 
about  270  hours)  is  characteristic.  The  “saturation  level”  is  different  for  crystals  of  dif¬ 
ferent  origin,  being  sometimes  as  low  as  50  percent  of  the  initial  transmission,  but  is  the 
same  for  sections  from  any  given  crystal. 
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Fig.  L  -  Typical  curves  of  L,iF  optical  transmission 
atl21.6mM  (the  strong  hydrogen  L.yman-a  line)  versus 
time  inair.  Crystals  were  cleaved  in  dry  argon  gas. 


Two  general  effects  suggest  themselves  as  possible  causes  of  the  observed  decrease 
in  transmission;  (a)  crystal  defects,  such  as  a  post-cleavage  surface  relaxation  of  the 
type  considered  by  Lennard-Jones  and  Dent  (2),  or  the  formation  of  a  space  charge  layer 
due  to  differing  energies  of  formation  of  positive  and  negative  ion  vacancies  as  calculated 
by  Lehovec  (3),  or  (b)  a  surface  layer  produced  by  adsorption  or  chemical  reaction. 

The  possibility  that  imperfections  were  responsible  w'as  examined  in  the  following 
manner:  Two  crystals,  cleaved  from  the  same  cube,  were  chemically  polished  (4)  after 
one  had  been  mechanically  polished  to  produce  in  it  a  dislocation  density  many  orders  of 
magnitude  higher  than  in  its  mate.  Both  crystals  showed  essentially  the  same  transmis¬ 
sion  and  decay  behavior.  The  surface  relaxation  mechanism  was  discarded  because  the 
decay  was  not  independent  of  environment,  as  is  shown  later.  Thus,  it  was  concluded  that 
mechanical  defects  were  probably  not  responsible. 

Since  mechanical  defects  did  not  seem  to  be  the  cause,  the  surface  layer  possibility 
was  investigated.  Moisture  is  a  common  enemy  of  the  alkali  halides,  and  LiF,  despite  its 
relatively  low  solubility,  is  apparently  no  exception.  This  is  demonstrated  in  Fig.  3  which 
shows  the  effect  on  two  crystal.s  which  were  exposed  to  room  air  for  250  hours  and  then 
transferred  to  the  wet  box.  The  high-humidity  exposure  clearly  accelerated  the  trans¬ 
mission  loss.  The  “saturation  level,”  however,  is  the  same  as  that  of  the  crystals  exposed 
to  room  air  (Fig.  2)  despite  continued  wet -box  exposure  to  600  hours. 

Additional  data  on  this  moisture  effect  are  shown  in  Fig.  4  for  two  thin  crystals  simul¬ 
taneously  cleaved  in  air.  One  was  measured  immediately  as  a  control  and,  as  curve  (a) 
shows,  was  quite  characteristic  in  its  initially  high  transmission,  decay  curve  (out  to  200 
hours)  and  approach  to  the  “saturation  level”  at  85  percent  of  the  initial  value.*  The  rest 
of  curve  (a)  will  be  discussed  in  a  later  section.  The  other  crystal  was  stored  55  days  in 
the  dry  box  before  being  measured.  The  results  of  subsequent  air  exposure  are  shown  in 
curve  (b).  The  step  discontinuities  are  due  to  vacuum  treatments,  discussed  later.  If  the 
crystals  were  truly  identical,  then  a  loss  of  about  6  percent  of  the  initial  transmission 
occurred  during  storage.  If  the  85  percent  saturation  level  characteristic  of  the  parent 


*A11  crystals  from  this  block,  unless  given  a  heat  treatment  or  chemical  polish,  eventually 
decayed  to  85  percent  of  the  initial  transmission. 
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Fig.  3  -  Optical  transmission  at  121.6  versus  time 
in  air.  Transfer  from  room  air  to  high-humidity  air 
at  250  hours  caused  more  rapid  decrease  to  "saturation 
level."  Crystals  were  cleaved  in  air. 


Fig.  4  -  Optical  transmission  at  121.6  m^  versus  time 
in  air.  (a)  at  200  hours,  a  20-hour  “good"  vacuum 
exposure  dropped  the  transmission  to  "saturation." 
Five  more  days  of  “good"  vacuum  had  no  effect.  At 
1042  hours,  3  hours  of  "good"  vacuum  at  500°C  caused 
further  decrease,  (b)  effects  of  “good"  vacuum  expo¬ 
sures  of  2  to  6  hours  are  shown  at  70,  135,  175,  and 
195  hours.  Crystals  were  cleaved  in  air. 
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crystal  block  is  assumed  for  curve  (b),  then  about  3  percent  loss  has  occurred.  Thus,  a 
relative  transmission  loss  of  less  than  6  percent  in  55  days  in  dry  air  as  compared  to  15 
percent  in  less  than  15  days  in  room  air  indicates  that  the  transmission  decay  is  very 
substantially  inhibited  by  dry-air  storage. 

The  further  step  of  replacing  air  by  an  inert  gas  was  also  examined.  Two  crystals 
were  cleaved  in  dry  argon  and  stored  there  for  46  days  before  the  data  in  Fig.  5  were 
taken.  The  control  curves  are  provided  by  Fig.  2  for  which  the  crystals  were  cleaved  at 
the  same  time  in  the  same  conditions.  Comparison  of  initial  values  shows  that  again  only 
a  small  decrease  of  around  3  percent  could  have  occurred  during  the  storage.”  Thus,  it 
concluded  that  moisture  is  the  principal  offender,  that  the  cleavage  atmosphere  has  no 
permanent  effect  since  both  crystals  decay  in  essentially  the  same  way  when  exposed  to 
room  air,  and  that  physical  changes,  such  as  surface  stress  relaxation  after  cleavage,  are 
not  responsible. 


Fig.  5  -  Optical  transmission  at  121.6  m^  versus  time 
in  air  for  two  LiF  crystals  cleaved  and  stored  in  dry 
argon  gas  for  46days  before  measurements.  Compari¬ 
son  with  control  curves  of  Fig.  2  indicates  less  than  5 
percent  loss  of  relative  transmission  during  storage. 


Since  this  is  a  surface  effect  primarily  dependent  on  the  atmosphere,  the  possibility 
of  an  adsorbed  layer  which  might  be  removed  by  reducing  the  ambient  pressure  was  tested. 
It  was  consistently  observed  that  the  transmission  of  samples  stored  in  the  monochromator 
(P  <  5^l)  did  not  change.  In  Fig.  2,  for  example,  the  first  7-2/3  hours  of  air  exposure 
actually  occurred  over  an  8 -day  period  in  segments  of  roughly  1  hour  each.  These  air 
exposures  were  separated  by  periods  of  24  hours  or  more  in  the  monochromator  during 
which  no  change  occurred.  Storage  in  a  “good”  vacuum  also  usually  had  no  effect.  How¬ 
ever,  prominent  exceptions  to  this  were  observed,  an  example  of  which  is  shown  in  Fig.  4, 
curve  (a),  at  200  hours.  At  this  point,  20  hours  in  a  “good”  vacuum  decreased  the  trans¬ 
mission  to  the  85  percent  “saturation  level.”  Continued  pumping  for  a  total  of  5  days 


*This  uncertainty  arises  from  the  fact  that  the  initial  transmissions  of  two  plates  cleaved 
from  the  same  block  can  differ  by  3  percent.  This  is  probably  caused  by  variations  in 
the  amount  of  plastic  bending  and  the  density  of  tear  markings. 
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caused  no  further  change.  This  saturation,  plus  the  lack  of  change  in  five  other  saturated 
crystals  simultaneously  exposed,  eliminated  oil  deposition  as  a  cause.  Smaller  decreases 
after  periods  of  2  to  6  hours  in  “good”  vacuum  are  also  shown  in  Fig.  4,  curve  (b),  and  in 
Fig.  5. 

In  order  to  check  this  vacuum-induced  decrease,  a  fresh  crystal,  cleaved  in  air,  was 
subjected  alternately  to  periods  of  0.5  to  7  hours  of  a  wet -box  exposure  and  of  3  to  40  hours 
in  a  “good”  vacuum.  The  accelerating  effect  of  the  wet  box  was  again  observed,  and  the 
“good”  vacuum  exposures  had  no  effect  until  after  43  hours  total  wet -box  exposure  (30  days 
after  cleavage).  At  this  point  the  vacuum-induced  decrease  was  observed,  dropping  from 
88  to  85  percent  in  2  hours  and  then  to  82  percent  in  100  hours.  Subsequent  wet -box  expo¬ 
sures  had  no  further  effect. 

An  explanation  of  this  apparently  anomalous  behavior  in  vacuum  may  lie  in  a  simple 
consideration  of  the  reaction  kinetics  involved.  Consider  the  hydrolysis  reaction  as  occur¬ 
ring  in  two  different  phases  distinguished  by  the  factors  determining  the  reaction  rate. 

In  the  first  phase  the  rate  is  governed  by  the  amount  of  water  present,  with  the  resction 
occurring  essentially  instantaneously  at  the  clean  crystal  surface.  During  this  phase, 
storage  in  dry  gas  or  vacuum  halts  the  reaction  simply  by  removing  the  water.  As  the 
reaction  proceeds,  a  film  of  reaction  products  builds  up  on  the  crystal  surface.  Incident 
water  molecules  must  now  diffuse  through  the  film  before  reaction  can  occur,  and  HF 
must  diffuse  out  again  to  maintain  stoichiometry.  The  vacuum-induced  decrease  occurs 
during  this  second,  diffusion -controlled  phase  as  a  result  of  more  rapid  removal  of  HF. 
Thus,  exposure  to  gaseous  HF  during  this  phase  might  be  expected  to  inhibit  the  decrease 
in  transmission  or  even  restore  some  of  the  loss. 

The  addition  of  heat  during  the  vacuum  exposure  had  the  effect  shown  in  the  last  sec¬ 
tion  of  curve  (a).  Fig.  4.  At  about  1042  hours  the  “saturated”  crystal  was  heated  (500”c 
for  3  hours)  in  “good”  vacuum.  Its  transmission  was  unaffected  by  the  treatment  itself 
but  dropped  rapidly  during  subsequent  air  exposure.  The  similarity  of  this  curve  to  those 
of  freshly  cleaved  crystals  (Fig.  2)  might  be  interpreted  as  being  due  to  similar  crystal 
surface,  except  for  its  substantially  lower  value. 

Several  attempts  were  made  to  observe  the  forward  scattering  of  light  from  surface 
roughening,  which  would  contribute  to  the  decreasing  transmission,  and  also  to  see  if  all 
regions  of  the  crystal  contribute  equally  to  the  diminution  in  transmission.  The  crystal 
was  mounted  in  the  spectrometer  light  beam  with  an  ultraviolet-sensitive  photographic  film 
close  behind  it.  A  number  of  crystals  were  thus  “photographed”  immediately  after  being 
cleaved  and  after  aging,  but  no  changes  due  to  scattering  were  detected.  The  most  strik¬ 
ing  feature  of  the  films  was  the  clear  appearance  of  the  surface  tear  markings  as  seen  in 
Fig.  6.  A  careful  study  of  these  markings  revealed  no  signs  of  broadening  with  age  which 
might  be  expected  to  occur  preferentially  at  surface  steps  if  the  surface  were  being 
attacked.  The  observed  effect  was  an  overall  diminution  in  the  transmission  of  the  crys¬ 
tal  with  no  increase  in  the  background  to  indicate  scattering  of  light. 

This  result  may  be  justified  from  a  knowledge  of  the  rate  of  decrease  in  transmission 
and  from  light  scattering  results.  Consider  the  simple  case  in  which  spheres  of  material 
of  a  differing  index  of  refraction  form  on  the  surface.  The  turbidity  coefficient  t  describes 
the  fractional  reduction  in  intensity  of  a  light  beam  in  passing  through  a  medium.  It  has 
been  shown  (5)  that  for  scattering  particles  much  smaller  than  the  wavelength  of  light,  the 
turbidity  is  proportional  to  the  sixth  power  of  the  particle  radius, 
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Fig.  6  -  A  transmission  photograph  of  a  cleaved 
LiF  crystal  made  with  1216A  light.  The  tear 
markings  on  both  crystal  surfaces  arc  clearly 
visible . 


where  R  =  radius  of  particle,  N  =  number  of  particles  per  unit  area,  and  t  =  time.  Thus 
in  order  that  the  observed  rate  of  change  of  transmission  be  attributable  to  light  scatter¬ 
ing,  it  would  be  necessary  that  dR/dt  decrease  at  least  as  fast  as  R"-’’,  a  result  which 
seems  unlikely.  For  example,  a  constant  rate  of  increase  in  sphere  volume  would  only 
make  dR/dt  vary  as  R"^. 

It  is  also  possible  that  the  decreased  transmission  is  caused  by  very  broad  or  over¬ 
lapping  absorption  bands  so  that  a  specific  band  is  not  seen.  Such  a  surface  film  can  be 
expected  to  show  interference  effects  with  the  first  maximum  or  minimum  occurring 
for  a  film  thickness  of  ■'/4n  where  is  the  wavelength  of  the  incident  light  (1216A),  and 
n  is  the  index  of  refraction  of  the  film.  Since  no  maximum  or  minimum  was  observed  in 
the  transmission  as  a  function  of  time,  an  upper  limit  could  be  placed  upon  the  film  thick¬ 
ness  if  the  surface  index  of  refraction  could  be  measured.  Such  a  change  in  index  of 
refraction  would  not  be  detected  by  the  usual  method  of  passing  light  through  a  prism  and 
measuring  the  deviation.  However,  the  relative  contributions  of  interference  and  absorp¬ 
tion  to  the  changing  transmission  could  be  determined  by  reflection  measurements,  since 
an  interference  loss  in  the  forward  direction  would  show  up  as  increased  reflection  in  the 
backward  direction. 

Since  permanent  prevention  and/or  removal  of  the  offending  surface  effect  seemed 
difficult  by  atmosphere  control,  the  possibility  of  modifying  the  crystal  itself  to  inhibit 
the  effect  was  briefly  examined.  Divalent  impurities  are  known  to  improve  the  cleavage 
properties  (6)  and  the  optical  transmission  (7)  of  other  alkali  halide  crystals,  both  of  which 
are  desirable  effects.  Hence,  a  crystal  of  LiF:Mg  was  grown*  by  using  random  pieces  of 
high-purity  LiF  from  Harshaw  Chemical  Company  as  raw  material,  with  increasing 
amounts  of  MgF^  added  at  regular  intervals.  The  resulting  crystal  cleaved  quite  easily 
but  showed  reduced  transmission  in  all  sections  and  this  approach  was  not  pressed  further. 

An  evaporated  layer  of  ]VIgF2  has  been  found  to  protect  metal  mirror  surfaces  (8)  and 
this  approach  was  also  tried  on  the  crystals.  A  layer  of  MgF^  about  40A  thick  was  evap¬ 
orated  onto  each  surface  of  several  crystals.  This  layer  caused  an  immediate  drop  in 
transmission  to  about  80  percent  of  the  precoating  value  (at  1216A).  The  transmission 
continued  to  decrease  with  air  exposure  at  a  more  rapid  rate  than  an  uncoated  crystal 


*The  crystal  was  grown  at  NRL  by  Mr.  William  Zimmerman,  formerly  of  the  Crystals 
Branch,  Solid  State  Division. 
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and  in  70  hours  had  reached  50  percent  of  the  precoating  value  and  was  still  going  down. 
These  results  indicate  that  MgF2  overcoating  is  detrimental  to  the  optical  transmission 
of  LiF  crystals  in  the  vacuum  ultraviolet. 

Thin  plates  of  LiF  can  be  successfully  cleaved  at  liquid  nitrogen  temperature,  but  this 
is  an  awkward  process  and  must  be  done  carefully  to  avoid  extensive  plastic  deformation 
during  the  thermal  cycle.  A  more  convenient  technique  is  to  harden  the  crystal  by  x- 
radiation.  A  marked  increase  in  ease  of  cleavage  was  obtained  by  a  1-hour  x-ray  exposure 
(225  kv,  10  ma,  tungsten  target).  The  effects  of  the  exposure  were  annealed  out  by  heat¬ 
ing  the  crystals  (500'’c,  dry  argon,  two  hours)  leaving  no  detectable  absorption  bands  in 
the  ultraviolet  from  3000A  to  1050A.  A  similar  result  has  been  recently  reported  by 
Nadeau  and  Johnston  (9)  who  found  that  an  exposure  of  1.6  x  10®  roentgens  at  room  tem¬ 
perature  made  it  possible  to  cleave  the  new  “soft”  LiF  and  that  the  radiation  damage, 
including  the  absorption  bands  at  445  and  242  m^,  could  be  removed  by  a  two-hour 
anneal  at  450°C.  Their  results  have  been  confirmed  in  this  investigation  and  extended  to 
show  that  no  radiation -induced  absorption  remains  in  the  ultraviolet  from  3000A  to  1050A 
after  annealing. 


SUMMARY 

It  has  been  found  that  a  general  decrease  in  transmission  below  1600A  occurs  when 
freshly  cleaved  LiF  crystals  are  exposed  to  room  air.  The  rate  and  amount  of  the  decrease 
vary  with  crystal  origin.  In  the  crystals  reported  herein  (all  from  the  same  block),  a 
total  decrease  at  1216A  of  about  15  percent  of  the  initial  transmission  value  was  essentially 
complete  in  about  10  days. 

The  principal  cause  of  the  decrease  is  a  surface  film  apparently  created  by  chemical 
reaction  with  moisture.  The  reaction  can  be  inhibited  by  early  storage  of  the  crystal  in 
dry  gas  or  in  vacuum.  After  some  loss  has  occurred  however,  it  is  accelerated  by  vacuum 
storage.  The  transmission  loss  can  be  partially  restored  by  chemical  removal  of  the 
surface  layer. 

It  was  found  that  cleavage  plates  suitable  for  use  as  optical  elements  in  the  vacuum 
ultraviolet  can  be  prepared  by  exposure  of  the  new  “soft”  LiF  to  ionizing  irradiation, 
followed  by  annealing  in  an  inert  gas  which  removes  the  radiation-induced  absorption. 
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SMALL  FATIGUE  CRACKS  AS  FRACTURE  ORIGINS 
IN  TESTS  OF  HIGH-STRENGTH  STEEL  SHEET 
[  Unclassified] 

J.  E.  Srawley 

The  Naval  Research  Laboratory  (Metallurgy  Division) 
is  investigating  the  crack  propagation  resistance  of  high- 
strength  steels  and  comparable  materials  for  use  in  highly 
stressed  structures.  Specimens  are  provided  with  cracks 
and  tested  in  tension  to  determine  the  dependence  of  tensile 
strength  upon  crack  dimensions  or  the  ability  to  endure 
repeated  yield  stress  loading,  and  the  results  are  inter¬ 
preted  in  terms  of  fracture  mechanics. 

Methods  were  developed  for  forming  small  semiellip¬ 
tical  surface  cracks  of  controlled  location  and  dimensions 
in  tensile  specimens  by  electrolytic  cracking  or  by  fatigue 
stressing  at  low  nominal  stress  levels.  The  exploitations 
of  these  experimental  methods  for  the  evaluation  of  steels 
proposed  for  construction  of  advanced  rocket  motor  casings 
will  complement  other  methods  of  evaluation  to  insure  that 
high-strength  steels  do  have  acceptable  crack  propagation 
resistance  at  yield  strength  levels  of  260,000  psi  and  higher. 

. FBM  Staff 


INTRODUCTION 

The  degree  of  ability  to  tolerate  small  cracks  is  an  important  factor  in  the  efficient 
utilization  of  high-stren^h  materials  stressed  in  tension.  The  ASTM  Special  Committee 
on  Fracture  Testing  of  High  Strength  Materials  has  the  responsibility  of  considering 
methods  of  evaluating  materials  in  this  respect  and  has  issued  two  reports  which  may  be 
consulted  for  general  guidance  (1,2).  Some  additional  information  is  to  be  found  in  a  more 
recent  publication  by  Irwin  and  Srawley  (3).  The  usual  type  of  specimen  employed  for 
crack  propagation  testing  of  high-strength  sheet  materials  is  provided  either  with  very 
sharp  edge  notches  or  with  a  central,  transverse  slot  with  very  sharp  ends.  It  has  been 
found  that  slots  with  fatigue  cracks  at  the  ends  are  most  effective.  Another  type  of  spec¬ 
imen  which  has  been  found  to  be  of  value  is  provided  with  a  small  surface  crack,  extend¬ 
ing  only  partway  through  the  thickness,  rather  than  with  a  slot  with  cracked  ends.  In  a 
previous  publication,  tests  of  specimens  in  which  the  cracks  were  formed  by  hydrogen 
embrittlement  cracking  were  discussed  (4).  It  was  shown  that  the  results  were  in  good 
agreement  with  a  relationship  proposed  by  Irwin  for  simielliptical  surface  cracks  (5)  and 
could  be  used  for  the  calculation  of  values  of  the  plane -strain  crack  toughness  Kj^. 

Fatigue  cracking  offers  an  alternative  method  of  introducing  small  surface  cracks 
into  sheet  specimens  which  has  certain  advantages  over  hydrogen  embrittlement  cracking 
or  stress  corrosion  cracking.  The  possible  effects  of  residual  hydrogen  or  corrodent  on 
the  results  of  the  subsequent  tests  are  eliminated,  and  also  the  fatigue  cracking  procedure 
is  more  convenient  to  perform  and  control.  However,  the  effective  sharpness  of  a  fatigue 
crack  may  depend  upon  the  conditions  of  its  formation.  Plastic  deformation  of  the  material 
in  the  vicinity  of  tlie  crack  front  occurs  during  fatigue  stressing,  the  extent  of  the  plasti¬ 
cally  deformed  region  being  greater  the  higher  the  ratio  of  the  maximum  nominal  stress 
to  the  yield  strength  of  the  material.  Because  of  the  many  complicating  factors,  there  is, 
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at  present,  no  way  of  predicting  just  how  the  fatigue  stressing  conditions  will  affect  the 
effective  sharpness  of  the  fatigue  crack.  It  is  therefore  necessary  to  make  experimental 
studies  which  will  provide  some  guidance  about  this. 


MATERIAL  AND  HEAT  TREATMENT 

The  specimens  were  obtained  from  a  single  sheet  of  consumable  electrode  vacuum- 
melted  steel  of  a  modified  AISl-SAE  H-ll-type  tool  steel  composition.  The  sheet  was 
nominally  1/16  inch  thick  and  contained  0.41%  C,  0.35%  Mn,  0.012%  P,  0.007%  S,  0.91%  Si,  4.81% 
Cr,  1.42%  Mo,  and  0.56%  V.  The  specimens  were  1.000  inch  wide  and  8.0  inches  long,  the 
lengths  being  in  the  rolling  direction  of  the  sheet. 

The  specimens  were  heat  treated  after  machining  to  size  and  prior  to  fatigue  crack¬ 
ing.  All  the  specimens  were  heat  treated  in  one  batch,  together  with  two  standard  sheet 
tensile  specimens.  The  heat  treatment  was  as  follows:  Austenitize  in  argon  at  1850±25°F 
for  30  minutes,  then  air  cool  individually;  temper  at  1000±10°F  for  2  hours  and  air  cool; 
finally,  the  temper  treatment  is  repeated  twice.  This  heat  treatment  resulted  in  an  ulti¬ 
mate  tensile  strength  of  263,000  psi,  a  yield  strength  of  213,000  psi,  and  an  elongation  on  2 -inch 
gage  length  of  6  percent.  From  previous  tests  of  this  material  it  was  known  that  the  crack 
toughness  would  be  low  in  this  condition.  This  relatively  brittle  condition  was  selected 
because  it  has  been  predicted  by  Irwin  that  sensitivity  to  crack  sharpness  in  a  crack  prop¬ 
agation  test  should  be  higher  the  more  brittle  the  material  and  the  higher  its  yield  strength  (3). 


FATIGUE  CRACKING 

A  schematic  diagram  of  the  fatigue  cracking  arrangement  is  shown  in  Fig.  1.  The 
specimen  is  clamped  as  a  cantilever  and  the  free  end  is  deflected  by  a  cam  rotating  at 
1750  rpm.  The  lower  support  tapers  to  a  point,  as  indicated  in  Fig.  1,  which  introduces  an 
indeterminate  stress  concentration  at  the  position  on  the  specimen  at  which  the  crack  is 
to  be  formed.  This  feature  allows  the  location  of  the  crack  to  be  controlled  and  lowers 
the  value  of  the  ratio  of  maximum  bending  moment  to  section  modulus  required  to  develop 
a  crack  in  a  practical  length  of  time.  Credit  for  the  idea  is  due  to  Yen  and  Pendleberry  (6). 
With  this  arrangement  it  is  possible  to  develop  fatigue  cracks  without  additional  stress 
concentration.  However,  it  was  found  that  cracks  would  be  initiated  much  more  rapidly 
if  small,  sharp  indentations  were  made  in  the  specimens,  as  indicated  in  Fig.  1.  These 
indentations  were  made  with  a  hardened  tool  steel  punch  prior  to  heat  treatment  of  the 
specimens. 

The  general  relationship  between  the  time  elapsed  before  carcks  were  observed  at 
the  ends  of  an  indentation  and  the  length  of  the  indentations  was  established  in  earlier 
work  (7)  and  is  illustrated  in  Fig.  2.  These  data  apply  to  a  different  steel  from  the  one 
used  in  the  present  work,  but  the  two  steels  had  almost  the  same  yield  strength.  On  the 
basis  of  Fig.  2  an  indentation  length  of  0.05  inch  was  selected  for  the  specimens  for  the 
present  work.  These  were  intended  to  have  fatigue  crack  lengths  of  about  0.10  inch  in 
every  case.  The  fracture  of  a  typical  specimen  is  shown  in  Fig.  3,  and  it  can  be  seen 
that  the  indentation  covers  only  a  small  part  of  the  area  of  the  fatigue  crack. 

For  the  arrangement  shown  in  Fig.  1  the  ratio  of  bending  moment  to  section  modulus 
is  given  by  6Eyh/{4A^  +  3LA),  where  E  is  Young’s  modulus,  y  is  the  deflection,  h  is  the 
specimen  thickness,  and  A  and  L  are  the  dimensions  indicated  on  Fig.  1.  This  is  equal  to 
230,000y  psi  for  the  constant  values  of  E,  h,  A,  and  L  that  applied  to  the  present  work. 
Three  sets  of  ten  specimens  each  were  fatigue  cracked  using  cams  which  applied  maximum 
deflections  of  0.263,  0.363,  and  0.563  inch,  respectively.  The  minimum  deflections  were 
zero  in  each  case.  Thus,  the  maximum  nominal  stresses  for  the  three  sets  were  60,500, 
83,500,  and  130,000  psi. 
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Fig 


INDENTATION  LENGTH  (IN.) 


Fig.  Z  -  Relationship  between  time  elapsed 
before  formation  of  cracks  at  the  ends  of 
the  indentations  and  indentation  length  for 
a  high-strength  steel.  These  data  are  for 
a  steel  different  from  the  one  used  in  the 
present  investigation,  but  the  two  steels  had 
very  similar  yield  strengths. 
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Fig.  3  -  Typical  fracture  of  a  &  1  x  1 /16-inch  specinaen.  Top: 

entire  fracture;  bottom:  enlarged  view  of  fatigue  crack.  Dimen¬ 
sions  shown  in  top  photo  are  1  inch  and  1/16  inch. 


During  fatigue  cracking  the  crack  lengths  were  measured  at  intervals  with  a  low- 
power  microscope  equipped  with  a  micrometer  eyepiece.  Figure  4  shows  the  dependence 
of  crack  length  on  the  number  of  cycles  for  each  of  the  three  sets  of  specimens.  This 
dependence  varied  among  the  ten  specimens  of  a  given  set,  and  Fig.  4  shows  the  range  of 
values  for  each  set.  It  is  clear  that  the  rate  of  crack  propagation  depends  strongly  on  the 
maximum  deflection,  the  average  crack  length  increase  per  thousand  cycles  being  0.0010, 
0.0026,  and  0.0142  inch  in  order  of  increasing  deflection. 


Fig.  4  -  Semilog  plot  of  overall  crack  length  vs  number  of 
cycles  at  which  the  specimen  was  fatigued  to  the  maximum 
nominal  stress  shown.  Each  band  represents  the  range  of 
values  for  the  ten  specimens  fatigued  at  that  stress  level. 
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RESULTS  AND  DISCUSSION 

The  fatigue-cracked  specimens  were  tested  in  tension  at  room  temperature  (75  to 
80“F)  in  the  laboratory  atmosphere.  The  lengths  and  depths  of  the  fatigue  cracks  were 
then  measured  with  a  traveling  microscope.  These  dimensions  and  the  gross  fracture 
strengths  of  the  specimens  are  given  in  Table  1.  Values  of  the  stress  intensity  factor  K 
were  calculated  from  these  data  in  the  manner  described  in  a  previous  publication  (4)  and 
are  also  given  in  Table  1. 

No  evidence  was  found  of  slow  crack  growth  occurring  at  a  load  appreciably  less  than 
the  maximum  load  in  any  of  the  tests.  A  crystal  pickup  was  attached  to  each  specimen 
during  the  tests,  but  no  noise  indicative  of  crack  extension  was  heard.  Furthermore, 
close  scrutiny  of  the  fracture  surfaces  failed  to  reveal  any  singularities  of  fracture  mark¬ 
ings  which  might  indicate  slow  crack  extension  (see  Fig.  3  for  instance).  Thus,  it  is  con¬ 
sidered  that  the  K  values  in  Table  1  correspond  closely  to  Kj^  plane-strain  crack  tough¬ 
ness  values  (2,3,4). 

Averages  and  standard  deviations  of  the  K  values  for  each  group  of  ten  specimens 
are  shown  in  Table  1.  The  averages  do  not  differ  significantly,  but  the  standard  deviation 
for  the  group  fatigued  at  60,500  psi  is  significantly  lower  than  those  for  the  groups  fatigued 
at  higher  stresses.  The  ratio  of  sample  variances  is  greater  than  would  occur  by  chance 
once  in  twenty  times. 

Examination  of  the  results  showed  that  the  relationship  between  the  K  values  and  the 
crack  dimensions  was  quite  random,  whereas  there  was  a  strong  correlation  between  the 
K  values  and  the  gross  fracture  stresses,  as  shown  in  Fig.  5.  There  are  three  factors 
which  might  contribute  to  the  variability  of  the  fracture  stress  and  K:  variability  of  mate¬ 
rial  toughness  from  specimen  to  specimen,  variability  among  the  fatigue  cracks,  and  pos¬ 
sible  variable  environmental  influences  during  the  tests.  The  third  factor  can  probably 
be  discounted  since  the  tests  were  all  conducted  under  the  same  conditions  and  lasted 
only  a  few  minutes.  It  is  likely  that  material  toughness  variability  accounts  for  some  of 
the  variation  in  the  results,  but  it  is  also  considered  likely  that  variability  among  the 
fatigue  cracks  is  the  major  factor,  on  account  of  the  significant  difference  in  stress- 
intensity-factor  standard  deviations. 

Additional  specimens  were  fatigued  and  then  sectioned  through  the  cracks  for  metal- 
lographic  examination.  No  obvious  differences  could  be  found  between  cracks  formed  at 
low  or  high  fatigue  stresses,  the  cracks  m  both  cases  being  very  fine  with  tips  beyond  the 
limit  of  resolution  of  the  optical  microscope.  The  question  of  possible  differences  in  the 
fine  structure  of  the  material  in  the  immediate  vicinity  of  the  crack  tips  had  to  remain 
unresolved  because  previous  experience  with  replication  electron  microscopy  of  the  mate¬ 
rial  had  shown  that  this  method  would  be  inadequate.  It  is  possible  that  information  might 
be  obtained  by  transmission  electron  microscopy,  but  the  problems  of  preparation  of  a 
suitable  “thin -film”  specimen  are  formidable. 

Apart  from  the  possibility  of  fine  structure  differences  due  to  the  difference  in  fatigue 
stress  levels,  there  is  also  a  strong  possibility  of  a  difference  in  residual  stresses.  The 
residual  stress  around  the  tip  of  a  fatigue  crack  might  vary  from  specimen  to  specimen 
and  might  be  expected  to  show  greater  variation  the  higher  the  fatigue  stress  level.  A 
test  of  this  hypothesis  could  be  made  by  fatigue  cracking  specimens  before  they  were  heat 
treated.  Presumably  any  residual  stresses  from  fatigue  cracking  would  be  relieved  dur¬ 
ing  tlie  heat  treatment,  and  any  due  to  heat  treatment  .should  be  less  variable. 
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Table  1 

Results  of  Tensile  Tests  on  Fatique -Cracked  Specimens 


Crack  Length 

Crack  Depth 

Gross  Fracture  Stress 

K 

(in.) 

(in.) 

(psi) 

(psi-in.^-^^) 

60,500  psi* 

0.1050 

0.0300 

154,000 

42,800 

0.1040 

0.0280 

159,700 

43,700 

0.1010 

0.2950 

159,500 

43,800 

0.1010 

0.2900 

165,000 

45,300 

0.1010 

0.2800 

168,500 

45,900 

0.1040 

0.0270 

161,000 

43,800 

0.1050 

0.0280 

165,000 

46,600 

0.1030 

0.0280 

161,700 

44,200 

42,800 

0.1020 

0.0270 

158,200 

0.1010 

0.0270 

170,500 

46,000 

Average  K:  44,550 

Standard  Deviation:  1,350 

83,500  psi* 

0.0995 

0.0265 

177,300 

48,000 

0.1015 

0.0275 

167,500 

45,600 

0.1015 

0.0272 

153,500 

41,400 

0.1024 

0.0277 

159,800 

43,500 

0.1004 

0.0268 

170,500 

46,200 

0.0994 

0.0268 

148,500 

39,500 

0.1006 

0.0279 

162,000 

44,100 

0.1000 

0.0270 

145,800 

38,800 

0.1020 

0.0271 

170,000 

46,400 

0.1000 

0.0274 

169,500 

46,200 

Average  K:  43,950 

Standard  Deviation:  2,960 

130,000  psi 

0.1034 

0.0300 

170,000 

47,400 

0.1080 

0.0318 

164,000 

46,600 

0.1038 

0.0298 

152,000 

41,900 

0.1006 

0.0295 

158,000 

43,300 

0.1002 

0.0291 

160,000 

44,000 

0.1012 

0.0291 

154,500 

42,100 

0.1022 

0.0299 

159,000 

43,800 

0.1022 

0.0288 

144,000 

39,200 

0.1022 

0.0289 

157,800 

43,000 

0.1020 

0.0289 

173,800 

47,800 

Average  K:  43,930 

Standard  Deviation:  2,580 

*Nominal  fatigue  stress  level. 
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GROSS  fracture  STRESS  (PSD 

Fig.  5  -  Plot  of  stress  intensity  factor  K  vs 
gross  fracture  stress  for  the  three  groups 
of  specimens  indicated  by  the  values  of 
fatigue  stress 


CONCLUSIONS 

Fatigue  stressing  in  bending  is  a  convenient  and  dependable  means  of  introducing 
small  surface  cracks  of  approximately  semielliptical  shape  into  sheet  or  plate  specimens. 
Plane-strain  crack  toughness  values  K  can  be  calculated  from  the  results  of  fracture 
tests  of  such  specimens  by  using  a  relationship  due  to  Irwin  (5),  The  average  K  values 
for  groups  of  ten  specimens  were  found  to  be  independent  of  the  fatigue  stressing  level, 
but  the  standard  deviation  was  found  to  be  significantly  lower  for  the  lowest  stress  level 
used  than  for  higher  stress  levels.  It  is  considered  that  this  difference  is  attributable  to 
factors  associated  with  the  fatigue  cracks  which  depend  upon  the  fatigue  stress  level.  Two 
possibilities  are  residual  stresses  and  the  fine  structure  of  the  material  in  the  vicinity  of 
the  crack  tip. 

The  results  indicate  that  the  nominal  fatigue  stress  level  at  which  a  crack  is  formed 
may  have  only  a  very  small  effect  on  the  result  of  a  subsequent  fracture  tou^ness  test, 
but  that  the  results  are  likely  to  be  most  consistent  when  cracks  are  formed  at  stress 
levels  which  are  as  low  as  can  conveniently  be  used. 
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APPLICATION  OF  FRACTURE  MECHANICS  TO  ADHESIVE  JOINTS* 

[Unclassified] 


E.  Riplingit  S.  Mostovoy,'}'  and  R.  Patrick} 

The  Naval  Research  Liaboratory  (Mechanics  Division), 
in  anticipation  that  the  performance  of  the  fleet  ballistic 
missile  will  be  materially  improved,  is  exploring,  on  a  long- 
range  basis,  the  possibilities  of  developing  new  materials 
cind  design  concepts  which  go  substantially  beyond  those 
available  now  for  the  fleet  ballistic  missile. 

In  order  to  assist  our  understanding  of  factors  which 
limit  the  strength  of  adhesive  bonds,  a  means  for  investigat¬ 
ing  these  bonds  within  the  analytical  framework  of  fracture 
mechanics  has  been  established. 

.....  FBM  Staff 


INTRODUCTION 

A  better  understanding  of  the  crack  resistance  of  adhesives  is  of  importance  in  fiber 
wound  rocket  chambers  for  two  reasons.  First,  there  is  the  obvious  need  for  information 
about  the  toughness  of  large  bond  areas  where  nozzles  and  bosses  are  attached  to  the  body 
of  the  chamber.  In  addition,  since  the  loads  externally  applied  to  the  chamber  are  dis¬ 
tributed  to  the  various  layers  of  windings,  and  indeed  to  individual  fibers  by  the  resin 
between  them,  the  adhesive  toughness  is  a  significant  factor  in  establishing  the  strength 
of  the  structure. 

The  present  project  is  concerned  with  evaluating  the  toughness  of  adhesives  by  using 
the  procedures  of  fracture  mechanics.  This  is  a  relatively  new  technique  which  has  met 
with  considerable  success  when  applied  to  homogeneous  systems  such  as  the  high-strength 
steels  used  for  pressure  vessels.  By  the  use  of  fracture  mechanics,  an  intrinsic  toughness 
of  the  material  is  measured  which  makes  it  possible  to  determine  the  load  carrying  capacity 
of  a  structure  in  the  presence  of  flaws.  Thus,  minimum  toughness  standards  can  be  estab¬ 
lished  for  the  structural  materials  used  in  a  specific  application. 

Fracture  mechanics  has  not  previously  been  applied  to  heterogeneous  systems  such 
as  adhesive  joints.  One  reason  for  this  is  that  the  difference  in  elastic  moduli  of  the 
systems’  components  makes  the  mathematical  analysis  required  for  calculating  stresses 
in  the  vicinity  of  a  crack  prohibitively  difficult.  To  avoid  the  need  for  carrying  out  these 
analytical  stress  analyses,  which  are  available  only  for  single  component  structures, 
testing  techniques  had  to  be  developed  that  would  permit  measurement  of  fracture  tough¬ 
ness  without  knowing  the  stresses  at  the  crack  tip. 

Adhesive  joints  also  differ  from  homogeneous  systems  in  that  cracking  may  occur  by 
a  number  of  different  modes.  While  a  metal  structure  will  almost  invariably  fail  due  to 
an  openii^  mode,  i.e.,  from  tensile  stresses  acting  normal  to  the  fracture  path,  adhesive 
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joints  may  also  fail  due  to  forward  or  edgewise  shear.  Thus,  the  experimental  testing 
procedure  must  enable  one  to  evaluate  toughness  for  each  of  these  modes  of  crack 
extension  as  well  as  the  opening  mode. 


EXPERIMENTAL  PROCEDURE 

One  method  that  G.  R.  Irwin  has  suggested  for  measuring  fracture  toughness  is  the 
use  of  a  “calibration”  bar  (1).  This  procedure  has  the  advantage  that  it  does  not  require 
any  information  about  the  distribution  of  stresses  at  the  crack  tip.  Instead,  it  measures 
the  total  strain  energy  lost  to  the  system  by  the  propagating  crack.  The  equation  for 
fracture  toughness  using  this  procedure  is 

d(l/M)  (j) 

2  da 

where  Q  =  strain -energy  release  rate,  p  =  applied  load  over  a  unit  of  specimen  width, 

(1/Mi  =  compliance  of  system,  and  a  =  crack  length.  When  P  is  the  critical  load  required 
for  the  crack  speed  dg/dt  ,then  from  Eq.  (1)  measures  the  resistance  or  critical  crack 
extension  force  for  the  particular  speed  da/d  t .  Furthermore,  Eq.  (1)  is  applicable  to 
any  mode  of  crack  extension.  If  the  applied  load  acts  normal  to  the  fracture  path  and 
compliance  is  measured  for  this  method  of  loading,  the  observed  fracture  toughness  is 
that  for  an  opening  mode  Q,p  .  If,  on  the  other  hand,  p^,  is  a  shear  force  applied  in  the 
direction  of  crack  extension  and  (I'M)  is  the  corresponding  compliance,  the  fracture 
toughness  measured  is  that  for  forward  shear  Evaluating  ;)  by  a  side  shearing 

mode  requires  the  application  of  a  shear  load  parallel  to  the  crack  front;  this  measure 
of  fracture  toughness  is  designated 


Method  for  Measuring 

To  measure  opening  mode  toughness,  a  splitting  force  was  applied  in  the  plane  of  the 
crack  path.  This  method,  rather  than  the  more  conventionally  used  procedure  of  remote 
loading,  permits  control  over  crack  growth.  Some  preliminary  experiments  indicated 
that  even  with  loading  in  the  crack  plane,  satisfactory  control  was  only  obtained  when  the 
compliance  of  the  adhesive  sample  was  large.  To  fulfil  these  requirements,  a  double 
cantilever  beam  specimen  of  the  type  shown  in  Fig.  1(a)  was  developed.  A  schematic 
diagram  of  the  loading  arrangement  and  a  description  of  the  testing  procedure  are  given 
in  Appendix  A.  The  test  is  carried  out  by  applying  a  load  through  fingers  that  fit  into 
hole  F.  Prior  to  collecting  data  on  adhesive  systems  a  master  curve  of  (l/M)  versus  a, 
and  from  this  a  second  curve  of  (i.-M)  versus  d(i/M)/da,  was  obtained  by  using  a  saw-cut 
solid  aluminum  sample.  With  the  latter  curve  it  was  only  necessary  to  collect  data  on 
versus  (i/M),  as  the  crack  was  extended,  to  calculate  from  Eq.  (1). 


Method  for  Measuring 

To  find  fracture  toughness  in  a  forward  slipping  mode,  the  specimen  in  Fig.  1(a)  was 
modified  as  shown  in  Fig.  1(b)  so  that  a  shear  load  could  be  applied  to  the  adhesive.  The 
test  piece  in  this  case  simply  consisted  of  a  tension  and  a  compression  member  on  either 
side  of  the  glue  line  so  that  d(i/M)/da  was  a  constant,  and  (Jjj  is  obtained  by  multiplying 
the  load  for  crack  extension  by  a  constant  which  is  independent  of  the  crack  length  a. 

Method  for  Measuring  and  Combined 

In  practice,  adhesive  systems  are  subjected  to  a  combination  of  tension  and  shear. 

To  evaluate  fracture  toughness  under  this  loading  condition,  the  specimen  is  tilted 
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on  the  loadit^  fingers  as  shown  in  Fig.  1(c).  The  tilt  angle  e  determined  the  ratio  of 
tensile  to  shear  load  developed  at  the  crack  tip. 


Method  for  Measuring 

Although  attempts  to  measure  Gmc  have  not  yet  been  made,  it  is  expected  that  this 
will  present  little  problem  using  the  specimen  shape  and  loading  method  shown  in  Fig.  1(d). 


Fig.  1  -  Methods  of  specimen  preparation  for 
measuring  (a)  ;Jj  (fracture  toughness  for  open¬ 
ing  mode),  (b)  (fracture  toughness  for  for¬ 
ward  shear),  (c)  combined  and  Qu  (9 
determines  the  ratio  of  tensile  to  shear  load 
developed  at  the  crack  tip),  and  (d)  ^III  (frac¬ 
ture  toughness  for  side  shearing  mode) 


Types  of  Adhesive  and  Adherends 

Only  a  single  adhesive  was  studied  to  date.  Because  at  the  initiation  of  the  study  it 
was  helpful  to  locate  the  crack  tip  by  observing  stress  fringes,  a  commercial  photostress 
material,  Budd  Photostress  Type  A,  was  used.  The  adhesive  has  an  elastic  modulus  of 
0.42  X  10®  psi,  a  tensile  strength  of  8850  psi  and  a  total  elongation  of  3.6  percent. 

Most  of  the  data  collected  was  on  an  adhesive  system  using  2024-T4  aluminum  as  the 
adherends.  The  method  of  preparing  samples  from  these  is  given  in  Appendix  B.  The 
glass  adherond  material  was  a  low-iron  crown  glass.  Its  characteristics  and  method  of 
preparation  of  samples  from  it  are  also  given  in  Appendix  B. 


TEST  RESULTS 

Effect  of  Joint  Thickness  on 

The  first  variable  to  be  examined  was  the  effect  of  joint  thickness  on  (Jj  using  alumi¬ 
num  adherends.  Fig.  2.  The  data  in  the  top  curve  were  collected  on  specimens  run  at  a 
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Fig.  2  -  Effect  of  bond  (joint)  thickness 
on  the  critical  crack  extension  force 
Qlc  for  two  strain  rates.  (Each  data 
point  is  an  average  of  approximately  20 
to  30  measurements.) 


dQ/dt  of  about  1  at  the  instant  the  crack  started  to  move.  This  corresponds  to  a  dis¬ 
placement  of  the  loading  fingers  of  approximately  0.01  inch/min. 

When  the  joint  was  very  thick,  in  the  range  of  40  to  90  mils,  the  epoxy  in  the  joint 
had  about  the  same  toughness  as  the  bulk  material.  Reducing  the  joint  thickness  caused 
a  drop  in  epoxy  toughness  to  a  minimum  at  about  20  mils.  A  further  lowering  of  the 
joint  thickness  caused  the  toughness  to  rise  abruptly  to  a  maximum  that  is  about  twice 
as  high  as  the  toughness  of  the  bulk  material.  Still  thinner  joints  caused  the  toughness 
to  drop  again.  The  toughness  peak  was  very  narrow,  approximately  one  mil  wide.  The 
exact  location  of  the  maximum,  as  well  as  its  width  and  height,  have  not  been  rigorously 
defined  as  yet.  Indeed,  on  the  basis  of  the  data  in  Fig.  3,  one  might  question  its  very 
existence.  However,  a  rise  to  about  the  same  level  or  higher  than  that  of  the  bulk  epoxy 
is  indicated  for  all  thickness  below  five  mils. 

The  difficulty  of  defining  the  exact  shape  of  the  Qj  curve  results  from  the  difficulty 
of  making  specimens  whose  joint  thickness  is  held  constant  along  the  entire  specimen 
length  and  width  and  the  further  problem  of  measuring  this  thickness  nondestructively. 

To  maintain  a  sample -making  px'ocedure  that  was  as  versatile  as  possible,  the  specimens 
were  made  without  the  benefit  of  a  jig,  and  as  a  consequence  there  were  significant  joint 
thickness  variations  in  all  the  samples.  With  thick  joints  this  presented  no  problems, 
but  in  the  thinner  joints,  it  obviously  was  serious. 


Fig-  3  -  Qjj,  measured  on  two  specimens 
with  tapered  joint  thickness  varying  from 
2  to  5  mils  along  the  specimen  length 
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The  occurrence  of  this  Qj  peak  was  veri!ieo-«ll>>7usiiig  samples  in  which  the  joint 
thickness  was  purposely  tapered.  Fracture  Ii)uf»iislii:a.e8s  results  obtained  on  two  samples 
whose  glue  line  nominally  varied  from  2  to  Sninills  along  the  specimen  length  are  shown 
in  Fig.  3.  Although  h^h  toughness  values  o\tr  iTair~estricted  range  of  thickness  were  found 
in  these,  the  thickness  that  produces  the  maximungicm  has  still  not  been  ascertained  even  in 
the  selected  specimen  shape.  Unfortunately eveinfen^when  the  location  of  the  peak  is  estab¬ 
lished,  it  is  expected  to  be  a  function  of  at  leiEtBstsp>ecinaen  width. 


Stability  Factor 


Although  (Jj^  is  an  adequate  measure  oflheaeet-iergy  req.ulred  for  crack  motion,  two 
joints  with  the  same  toughness  may  have  differeaejt  fracturing  behaviors  depending  on  the 
geometry  of  the  system.  This  results  becauKtsIbe?  strain  energy  in  the  system  is  stored 
in  the  adherends.  As  energy  is  consumed  bjthil-tea_dvancing  crack,  a  “stiff"  system  will 
consume  only  a  small  amount  of  its  total  eneigyKjr-'hile  a  “soft"  system  will  release  a 
significant  portion  of  its  energy.  Crack  groiththin  the  latter  then  will  be  far  more  stable 
than  in  the  former.  In  order  to  quantitativeljm  rela_te  joint  configurations  with  ability  to 
arrest  crack  propagation,  G.  Irwin  (2)  definrf  “stability  factor"  for  a  stable  crack 
extension  tendency  as  the  fractional  decreaseiri  iii5  per  unit  of  crack  extension  under  fixed 
grip  conditions,  i.e.,  as 


He  suggested  that  when  the  stability  factor  vss  i  s less  than  about  0.5  in.'‘  propagation  of 
cracks  takes  place  uncontrollably.  Using  thisd  ddSinition,  the  stability  factor  has  been 
calculated  for  specimens  with  various  heighisw  ho®’  adherends  (Fig.  4).  Note  that  the 
specimens  become  increasingly  stable  as  h  isd  te  xeased.  This  chart  would  indicate  that 
a  3-inch-high  specimen  has  about  the  same  stab  ahilfi-ty  with  a  6-1/2-inch  crack  as  a  1/2-inch 
sample  has  with  an  8-inch  crack.  Although  weoexp^erimental  data  are  available  to  support 
these  curves,  the  behavior  is  generally  consistdstwt  with  the  authors’  experience  in  Hpveloping 
the  test  sample  shape. 


The  crack  length  at  which  control  could » *-olotnger  be  exercised  over  crack  growth  was 
noted  for  the  specimens  used  to  get  the  effect oJ  ofjaint  thickness  on  Qj .  The  results  are 
shown  in  Fig.  5.  This  would  indicate  that  thepi!  pi«p>osec3  expression  for  stability  factor 
might  better  be  written  as 


1  L5 

h  (Jda  'i-a 


(3) 


Expression  (3)  indicates  a  hyperbolic  relatiissBi-sliip  between  joint  thickness  and  stability 
factor  as  defined  in  Eq,  (2).  Since  the  bulk  epsfnxjt^ecame  unstable  with  a  five -inch -long 
crack,  this  seems  a  reasonable  expression  fcr  trlbi,  ck  joints .  The  hyperbolic  relationship 
is  not  evident  in  thin  joints  because  the  crack  W  length  at  instability  exceeds  the  specimen 
length. 


Effect  of  Strain  Rate 

A  very  slow  strain  rate  was  used  to  colleclisclttxe  data  discussed  above;  however,  a  few 
tests  have  also  been  run  at  a  faster  rate.  Ttie*selfi^ctof  the  increased  strain  rate  on  the 
fracture  toughness  of  bulk  epoxy  is  shown  bjthWbepiJroperty  loss  in  the  solid  material  on  the 
extreme  right  in  Fig.  2.  Again  with  thick  joiat^*ts,t:lie  epoxy  when  used  as  an  adhesive  had 
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CftACK  LENGTH  0  (INCHES) 

Fig.  4  -  Stability  factor  i/(J  d(J/da  as  a  function  of  crack  length 
for  specimens  of  varying  height  li 


BONO  THICKNESS  (MILS) 

Fig.  5  -  Dependence  of  crack  instability 
on  joint  thickness 
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the  same  toughness  as  the  bulk  material,  and  the  strain  rate  effect  appeared  to  be  inde¬ 
pendent  of  joint  thickness  down  to  about  20  mils.  Sufficient  tests  have  not  yet  been  run 
to  describe  the  effect  of  strain  rate  on  thinner  joints. 


Effect  of  Joint  Thickness  on 

In  a  limited  number  of  tests,  was  measured  on  samples  of  varying  joint  thickness. 
The  data,  shown  by  circles  in  Fig.  6,  are  not  sufficiently  complete  to  define  the  effect  of 
joint  thickness  on  Qn  ,  but  in  all  cases  the  joint  is  at  least  an  order  of  magnitude  tougher 
in  shear  than  in  tension.  The  slow  strain  curve  of  Fig.  2  is  added  to  this  chart  for 
comparison. 


Fig.  6  -  Effect  of  joint  thickness  on  '*IIc 
(circles).  The  values  of  from  Fig.  2 
are  shown  for  comparison. 


Effect  of  Combined  Tension  and  Shear  on  (Jj,, 

A  few  preliminary  measurements  of  (Ji  with  a  shear  load  superimposed  onto  the  tensile 
load  were  also  made.  In  all  cases  the  angle  of  tilt  O,  shown  in  Fig.  1,  was  45  degrees  so 
that  identical  shear  and  tensile  loads  were  developed  at  the  crack  tip.  Surprisingly,  the 
addition  of  the  shear  component  onto  the  tensile  load  results  in  an  appreciable  increase  in 
opening  mode  joint  toughness,  as  shown  in  Table  1. 


Table  1 

Effect  of  a  Superimposed  Shear  Load  on 


Joint 

Thickness 

(mils) 

Ratio  of 
Shear  Load 
to  Tensile  Load 

Strain  Rate 

4ic 

B 

20 

0 

1 

1.3  (avg) 

20 

1 

1 

4.6 

10 

0 

25 

0.8  (avg) 

10 

25 

1.7 
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Effect  of  Adberend 

The  results  described  above  were  all  obtained  on  an  aluminum-epoxy-aluminum  system. 
Since  aluminum  and  glass  have  the  same  elastic  modulus,  bonds  made  with  aluminum  and 
glass  adherends  would  be  expected  to  have  the  same  toi^hness  as  long  as  identical  adhesive 
and  joint  configuration  are  used.  On  testing  glass-epoxy-glass,  however,  values  were 
obtained  that  were  far  lower  than  those  found  with  aluminum  adherends  (Table  2).  These 
lower  values  of  toughness  could  be  traced  to  poor  adhesion  of  the  selected  epoxy  to  the 
glass  as  compared  with  its  adhesion  to  aluminum.  When  a  water  drop  was  placed  on 
opposing  fracture  surfaces  of  the  aluminum  samples,  the  two  contact  angles  were  the 
same  (approximately  60  degrees).  On  the  other  hand,  the  contact  angles  measured  on 
the  glass  adherend  fracture  surfaces  were  quite  different  (35  to  40  degrees  on  one  surface 
and  60  degrees  on  the  other).  Initially,  organic  cleaning  was  used  on  the  glass.  Acid 
cleaning  was  found  to  lower  the  toughness  while  flame  cleaning  improved  it.  Even  with 
the  latter,  the  toughness  of  the  bonds  using  glass  adherends  was  only  about  one -third 
that  of  the  joints  made  with  aluminum. 

The  cleaning  procedures  used  on  both  the  aluminum  and  glass  samples  is  given  in 
Appendix  B.  In  addition  to  the  difference  in  cleaning  method  used  on  the  two  adherends, 
the  bonded  aluminum  surfacp  was  considerably  rougher  than  the  glass,  as  described  in 
Appendix  B. 


Table  2 

Comparison  of  Toughness  Obtained  on 
Aluminum -Epoxy- Aluminum  and 
Glass-Epoxy-Glass  Systems  (Oij.  =  25) 


Adherend 

Joint 

Thickness 

(mils) 

Method  of 
Surface  Cleaning 

^Ic 

Aluminum 

10 

acid 

1.4 

5 

acid 

2.5 

Glass 

10 

organic 

0.33 

10 

acid 

0.24 

10 

flame 

0.47 

5 

organic 

0.23 

Fracture  Appearance 

The  appearance  of  the  fracture  surface  depended  on  the  type  of  adherend,  the  joint 
thickness,  and  the  stress  state  under  which  the  specimen  was  fractured.  Thick  bonds, 
using  aluminum  adherends  cracked  by  an  opening  mode,  invariably  fractured  in  the 
middle  of  the  epoxy  layer.  The  fracture  surface  was  smooth,  only  interrupted  by  a 
number  of  finger-nail  markings  along  the  specimen  length,  presumably  at  crack  arrest 
position.  As  the  bonds  were  made  thinner,  the  fracture  surface  generally  was  rougher. 
The  glass  samples,  on  the  other  hand,  always  broke  near  one  interface,  with  a  completely 
smooth  surface  having  no  interruptions. 

A  distinctive  crack  pattern  was  formed  in  the  aluminum  samples  broken  by  forward 
shear.  The  fracture  was  close  to  one  interface  in  these  samples  and  a  series  of  45-degree 


NAVAL  RESEARCH  LABORATORY 


47 


Fig.  7  -  Typical  appearance  of  a  (Jjj  specimen 


cracks  formed  in  the  epoxy,  indicating  that  the  failure  was  the  result  of  a  normal  tensile 
stress.  The  typical  appearance  of  a  Qjj  specimen  is  shown  in  Fig.  7. 
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APPENDIX  A 


METHOD  FOR  COLLECTING  AND  DATA 


A  schematic  drawing  of  the  device  used  for  collecting  the  data  required  for  (Jic  calcu¬ 
lations  is  shown  in  Fig.  Al.  A  linear  variable  differential  transformer  is  mounted  on  the 
specimen  near  the  loading  holes,  after  which  the  specimen  is  placed  on  the  calibrated 
loading  fingers  and  the  transducer  is  calibrated  with  shims  for  accuracy  of  vertical  dis¬ 
placement.  Using  a  rectified  variable  dc  power  supply,  the  drive  motor  is  activated  by 
applying  a  voltage  corresponding  to  a  given  strain  rate  that  has  been  previously  deter¬ 
mined.  The  outputs  of  both  the  strain  gage  load  cell  and  the  transducer  are  recorded 
simultaneously  on  an  x-y  recorder  with  its  ranges  set  for  best  determination  of  slope 
and  critical  load  for  a  given  crack  length. 

A  microscope  and  polarized  light  source  are  also  mounted  such  that  the  crack  can 
be  observed  and  its  length  measured.  Because  it  is  difficult  to  accurately  determine  the 
exact  location  of  the  crack  tip,  the  observed  crack  lengths  are  not  used  in  calculating 
Instead,  a  calculated  crack  length  determined  from  the  calibrated  compliance  of  the  system 
is  used. 

If  the  crack  is  controllable  and  the  loading  rate  sufficiently  slow,  the  crack  is  allowed 
to  run  for  about  1/4  inch  between  measurements.  Once  the  crack  has  run  a  short  distance, 
the  motor  is  electrically  reversed  until  the  specimen  is  unloaded.  This  procedure  is 
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Fig.  Al  -  Schematic  of  tensile  loading  arrangement 
with  enlarged,  detailed  view  of  the  split  circular 
membrane  (crack-opening  device) 
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repeated  approximately  20  to  30  times  for  a  single  specimen  until  the  crack  becomes 
unstable.  The  critical  loads  and  slopes,  as  well  as  crack  length  at  instability,  are  logged 
for  calculations. 

An  identical  procedure  is  used  for  collecting  data  with  a  superimposed  shear  load. 

When  the  crack  is  propagated  by  forward  shear,  the  test  specimen  consists  simply  of 
one  tension  and  one  compression  member.  Thus  the  calculated  value  of  d(i/M)/da  is  a 
constant  and  the  evaluation  of  only  i  aires  measurement  of  the  load  at  which  the 
crack  runs.  Because  the  adhesive  is  so  tough  in  shear  as  compared  with  tension,  it  was 
necessary  to  use  a  hydraulic  tensile  machine  in  a  range  capable  of  5000  pounds  for  col¬ 
lecting  these  data.  The  load  is  applied  by  gripping  one-half  of  the  specimen  in  the  V-grips 
of  the  machine  and  applying  a  compressive  load  to  the  other  half  of  the  specimen  by  means 
of  a  yoke. 


APPENDIX  B 


CLEANING  PROCEDURE  FOR  ALUMINUM  AND  GLASS  SAMPLES 


Method  of  Preparing  Aluminum  Samples 

The  aluminum  plates  used  as  adherends  are  milled  and  ground  on  the  12  x  i/4-inch 
face  with  an  open  wheel,  after  which  the  pieces  are  cleaned  as  follows: 

1.  Degrease  in  trichloroethylene 

2.  Vapor  degrease  with  trichloroethylene 

3.  Treat  the  metal  part  by  immersion  for  10-12  minutes  at  150°-170‘’F  in  an 
etching  solution  consisting  of: 


a.  Sodium  dicljroraate,  Na2Cr207-2H20  .  200  pbw 

b.  Sulfuric  Acid  (sp.  gr.  1.84)  .  700  pbw 

c.  Distilled  Water  .  1700  pbw 


4.  Rinse  thoroughly  in  cold  tap  water 

5.  Final  rinse  in  distilled  water 

6.  Dry  in  forced-draft  oven;  the  metal  temperature  should  not  exceed  140°-150®F 

Two  of  these  cleaned  pieces  are  placed  on  a  leveling  plate,  and  the  surfaces  to  be 
bonded  are  separated  by  a  shim  of  the  required  joint  thickness  placed  at  the  two  ends  of 
the  gap.  A  plug  is  placed  into  the  loading  hole  and  Teflon  tape  is  pressed  over  the  length 
of  the  gap  between  the  plates  to  form  a  dam  for  accepting  the  epoxy.  The  assembly  is  then 
placed  on  a  hot  plate  at  150"?  and  held  for  45  minutes  to  attain  equilibrium. 

The  epoxide  components  are  individually  warmed  to  105 "F  and,  just  prior  to  pouring, 
the  two  components  are  thoroughly  mixed  and  the  temperature  of  the  mixture  brought  to 
105"F.  The  bond  is  filled  by  pouring  through  a  polyethylene  funnel  having  a  l/32-inch 
orifice.  Pouring  started  at  the  plugged  loading  hole  and  the  funnel  is  drawn  along  the 
length  of  the  sample  at  a  speed  such  that  the  adhesive  overflows  the  bond  and  maintains 
a  bank  of  adhesive  in  front  of  the  tip.  This  procedure  avoids  air  entrapment  and  results 
in  bubble -free  samples. 


Method  of  Preparing  Glass  Samples 

The  mating  surfaces  of  the  glass  adherends  were  ground  to  seven  wavelengths  of  light. 
Table  B1  summarizes  the  properties  of  the  glass. 

Table  B1 

Properties  of  Low-Iron  Crown  Glass 


=  1.523;  V  =  58.8  (dispersion) 

Coeff.  of  Expansion  =  9.4  x  10*®/°C  (25°-300°C  range) 
Softening  Point  =  1343°C 

Color  Free  Transmission  =  A/88%  illuminant  at  0.5  in.  thick 
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The  three  alternate  cleaning  procedures  used  were  (a)  organic  cleaning,  (b)  acid 
cleaning,  and  (c)  flame  cleaning.  These  procedures  are  summarized  in  Table  B2. 


Table  B2 

Methods  for  Preparing  Glass 

I.  a.  Wash  thoroughly  with  detergent  solution 

b.  Dry  with  hot  air  gun 

c.  Wash  with  triply  distilled  benzene 

II.  a.  Wash  thoroughly  with  detergent  solution 

b.  Dry  and  place  in  mixture  of  sulfuric  acid 
(90)/nitric  acid  (10)  for  30  minutes 

c.  Rinse  thoroughly  with  distilled  water 

d.  Dry  with  hot  air  gun 

in.  a.  Wash  thoroughly  with  detergent  solution 

b.  Dry  with  hot  air  gun 

c .  Wipe  off  with  benzene 

d.  Gently  flame  with  an  oxidizing  flame 


The  samples  are  prepared  immediately  prior  to  assembling  in  order  that  the 
surfaces  will  not  become  contaminated.  Within  30  minutes  of  cleaning,  the  glass-epoxy¬ 
glass  specimens  were  poured  and  cured. 
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A  NEW  SUBMARINE  DETECTION  SYSTEM 
[Unclassified  Title] 

A  NEW  SYSTEM  FOR  DETECTION  OF  RADIOACTIVITY 
INDUCED  IN  NUCLEAR  SUBMARINE  WAKES 
[Secret  Title] 


G.  A,  Burginyon,  J.  W.  Butler,  and  W.  H.  Jones 

The  Naval  Research  Laboratory  (Nucleonics  Division) 
is  presently  engaged  in  a  programto  increase  the  knowledge 
of  the  environment  in  which  a  submarine  operates  and  to 
develop  a  method  of  obtaining  information  pertaining  to  the 
influence  of  the  submarine  on  this  environment. 

Previous  efforts  toward  the  detection  of  the  radioac¬ 
tivity  left  in  the  wake  of  a  nuclear -powered  submarine  have 
involved  the  use  of  simple  ungated  Geiger  counters  or 
ungated  scintillation  detectors.  The  system  described 
herein  utilizes  a  time -coincidence  scheme  to  reduce  normal 
background  with  respect  to  the  counting  rate  expected  from 
the  residual  radioactive  isotopes  sodium-24  and  chlorine-38 
from  the  wake  of  the  nuclear -powered  submarine. 

. FBM  Staff 


INTRODUCTION 

The  reactor  in  a  nuclear  submarine  is  a  prolific  source  of  neutrons,  some  of  which 
penetrate  the  hull  and  enter  the  surrounding  water.  The  great  majority  of  these  neutrons 
are  captured  by  elements  of  the  sea  water  which  lead  to  no,  or  very  weak,  radioactivity. 
However,  a  small  fraction  of  the  neutrons  is  captured  by  sodium  nuclei  producing  Na*^ 
which  has  a  half-life  ox  about  15  hours.  Another  small  fraction  of  the  neutrons  is  captured 
by  some  of  the  chlorine  nuclei  producing  Cl^,  which  has  a  half-life  of  about  37  minutes. 
The  research  to  be  described  involves  the  detection,  by  means  of  radiation -sensitive 
instruments,  of  the  Na^'’  and  Cl^  isotopes  in  the  wake  of  a  nuclear -powered  submarine. 

The  concept  of  detecting  the  radioactive  trail  of  a  nuclear  submarine  is  not  new.  The 
original  specifications  for  nuclear  submarines  included  a  maximum  permissible  neutron 
leakage,  calculated  to  make  the  radioactivity  of  sodium  in  the  wake  small  compared  with 
the  natural  radiations  in  the  sea,  mainly  from  so  that  it  would  be  essentially  impos¬ 
sible  to  detect  the  wake.  A  number  of  subsequent  independent  studies  have  also  shown 
that  the  detection  of  the  radioactive  trail  is  indeed  impossible  with  conventional  techniques. 

The  system  to  be  described  differs  from  conventional  systems  by  taking  advantage 
of  a  particular  characteristic  of  Na^**  and  Cl^  which  enables  one  to  differentiate  them 
from  the  much  more  intense  natural  radiation  in  the  sea  water.  The  characteristic 
referred  to  is  that  both  Na^'*  and  Cl^”  emit  two  gamma  rays  simultaneously  when  one 
such  nucleus  decays,  whereas  the  radioactive  atoms  naturally  occurring  in  sea  water, 
for  example  K"*® ,  emit  only  one  gamma  ray  for  each  decay.  Thus,  by  having  two  inde¬ 
pendent  detectors  of  gamma  rays  and  by  recording  only  those  events  when  both  detectors 
count  at  the  same  instant,  one  can  virtually  eliminate  the  background  due  to  and  most 
other  natural  radioactive  atoms  in  sea  water.  The  detection  of  a  radioactive  trail,  which 
had  formerly  appeared  impossible,  now  shows  promise. 


SECRET 


53 


54 


NAVAL  RESEARCH  LABORATORY 


SECRET 


APPARATUS  AND  PROCEDURE 

A  water  sample  is  removed  from  the  wake,  and  the  salt  is  concentrated  by  evaporation 
and  then  placed  in  a  thin  rubber  container  between  two  detectors,  which  are  sodium  iodide 
(thallium  activated)  crystals  5  inches  in  diameter  and  6  inches  long.  A  block  diagram  of 
the  two-dimensional  coincidence  analyaer  is  shown  in  Fig.  1.  Except  for  the  oscilloscope 
driver  circuits  and  the  scaler,  the  system  is  completely  transistorized. 

There  are  two  crystals  with  identical  parallel  circuits  detecting  gamma  rays  emitted 
by  the  radioactive  sample.  Let  us  suppose  that  a  single  gamma  ray  from  the  sample 
interacts  in  crystal  y.  The  3 -inch  multiplier  phototube  (type  6363)  optically  coupled 
to  the  crystal  will  send  a  small  pulse  to  the  gamma-ray  amplifier.  Here  the  pulse  is 
amplified  and  then  fed  both  to  a  Schmitt  trigger  circuit  and  to  a  delay  line.  The  delayed 
pulse  is  further  amplified  in  the  oscilloscope  driver  circuit  to  a  maximum  of  180  volts 
and  applied  directly  to  the  vertical  deflection  plates  of  the  cathode-ray  tube  (3RP11A). 

The  other  amplifier  output  connects  to  a  Schmitt  trigger  circuit  which  emits  a  constant 
size  and  shape  pulse  for  input  pulses  above  a  certain  fixed  amplitude.  This  pulse  is  now 
fed  to  a  fast  coincidence  circuit,  consisting  of  biased  diodes,  which  will  produce  no  output 
pulse  when  only  a  single  input  pulse  is  applied.  Since  the  output  from  the  fast  coincidence 
circuit  is  used  to  turn  on  the  beam  of  the  cathode-ray  tube  (which  is  normally  turned  off), 
no  trace  can  occur  on  the  oscilloscope  for  a  single  event  in  either  crystal.  If,  however, 
two  pulses,  one  from  each  of  the  two  crystals,  arrive  at  the  fast  coincidence  circuit  within 
0.1  microsecond  of  each  other,  the  circuit  will  be  activated  and  a  pulse  will  be  sent  to  the 


Fig.  I  -  Block  diagram  of  the 
coincidence -anticoincidence 
electronic  detector  system 
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oscilloscope  driver  circuit.  Here  the  pulse  is  effectively  amplified  and  shaped  and  then 
used  to  turn  on  the  beam  of  the  cathode -ray  tube.  Thus  if  two  gamma  rays  interact 
simultaneously  in  the  two  crystals,  one  gamma  ray  in  crystal  y  and  one  gamma  ray  in 
crystal  x,  then  simultaneous  pulses  will  occur  in  both  halves  of  the  circuit.  These  pulses 
will  activate  the  fast  coincidence  circuit,  which  turns  on  the  cathode-ray  beam,  and  they 
will  also  appear  on  the  deflection  plates  of  the  cathode-ray  tube,  thus  producing  a  trace 
on  the  oscilloscope.  The  voltage  on  the  vertical  deflection  plates  of  the  oscilloscope  is 
proportional  to  the  energy  deposited  in  crystal  y  and  the  voltage  on  the  horizontal  deflec¬ 
tion  plates  of  the  oscilloscope  is  proportional  to  the  energy  deposited  in  crystal  x.  It 
is  thus  possible  to  calibrate  the  ordinate  and  abscissa  of  the  oscilloscope  in  Mev.  There¬ 
fore,  for  gamma  rays  which  are  detected  simultaneously,  the  oscilloscope  will  show  a 
trace  which  terminates  at  the  coordinate  (E„  Ey),  where  E„  corresponds  to  the  energy 
deposited  in  crystal  x  and  Ey  corresponds  to  the  energy  deposited  in  crystal  y. 

The  double  photopeak  interaction  of  the  Na^*  isotope  (which  results  when  all  the 
energy  of  each  gamma  ray  is  deposited  in  its  respective  crystal)  will  produce  a  compact 
cluster  of  traces  terminating  in  a  small  area  at  coordinate  (E^,  Ej)  and  a  similar  cluster 
at  coordinate  (Ej,  Ej),  where  Ej  and  Ej  are  the  energies  of  the  Na^^  gamma  rays  (1). 

A  cluster  of  traces  will  also  appear  at  coordinates  (E3,  E^)  and  (E^,  Eg)  caused  by  the 
Cl^  gamma  rays.  (See  Fig.  2  for  a  simplified  diagram  of  the  decay  schemes  of  Na^^  and 
Cl'^.)  The  reason  for  the  two  concentrations  of  traces  fur  each  radioactive  isotope  is 
that  either  gamma  ray  can  enter  either  crystal,  resulting  in  symmetry  about  a  45 -degree 
line  from  the  origin  of  the  trace  provided  that  the  gains  in  both  halves  of  the  circuit  are 
equal. 


Fig.  2  -  A  simplified  diagram  of  the 
decay  schemes  of  the  Na  and  Cl-*® 
isotopes.  Only  those  features  in  the 
decay  scheme  pertinent  to  the  pres¬ 
ent  report  are  included. 


(S'*  No” 


In  addition  to  double  photopeak  interactions  there  are  Compton-photopeak  interactions 
and  Compton -Compton  interactions  in  the  two  crystals.  A  Compton  interaction  usually 
leads  to  the  situation  where  only  a  part  of  the  gamma-ray  energy  is  deposited  in  the 
crystal,  the  rest  of  the  energy  leaving  the  crystal  in  the  form  of  a  lower  energy  gamma - 
ray  photon.  (A  schematic  diagram  of  a  typical  display  on  the  oscilloscope  face  is  shown 
in  Fig.  3.  In  order  to  avoid  confusion  only  double  photopeak  traces  are  shown  in  the 
display.) 
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Fig.  3  -  Schematic  representation  of 
the  display  produced  on  the  oscillo¬ 
scope  face  by  the  Na*^  and  Cl^*  iso¬ 
topes.  To  avoid  confusion  only  those 
traces  resulting  from  a  double  pho¬ 
topeak  interaction  are  shov/n.  The 
dotted  lines  enclose  the  area  (Area  II) 
in  which  most  traces  caused  by  Na** 
and  Cl^®  terminate.  About  90  per¬ 
cent  of  the  background  traces  termi¬ 
nate  in  Area  I.  Occasional  high- 
energy  background  traces  will  sweep 
through  Area  II  and  terminate  in 
Area  III.  The  coordinate  (E^,,  E^) 
for  the  two  possible  locations  of  the 
trace  cluster  for  each  isotope  are 
also  indicated. 


The  face  of  the  oscilloscope  is  divided  by  dotted  lines  into  three  regions  labeled 
Area  I,  11,  and  III  in  Fig.  3.  About  90  percent  of  the  background  traces  terminate  in 
Area  I,  most  Na^'’  and  Cl^*  double  photopeak  and  Compton -photopeak  traces  terminate 
in  Area  II,  and  high-energy  background  traces  due  mostly  to  cosmic -ray  muons  termi¬ 
nate  in  Area  III. 

Two  three-quarter -inch  phototubes  (6365)  are  optically  coupled  to  Area  II  by  a  lucite 
light  pipe.  Two  additional  phototubes  are  optically  coupled  to  Area  in  in  a  similar  manner. 

If  a  trace  terminates  in  Area  ll,  the  phototubes  watching  this  area  produce  an  output  pulse. 
This  pulse  is  amplified  and  fed  to  a  Schmitt  trigger  circuit  where  it  is  shaped  and  sent 
through  an  anticoincidence  circuit.  The  pulse  is  now  fed  to  a  cable  driving  circuit  and 
finally  to  a  scaler  where  a  count  is  recorded.  If,  however,  a  trace  sweeps  through  Area  II 
and  terminates  in  Area  in,  a  rejection  pulse  will  be  produced  by  the  phototubes  watching 
Area  III,  in  addition  to  the  pulse  from  Area  II  which  was  described  above.  After  ampli¬ 
fication  followed  by  shaping  in  a  Schmitt  trigger,  this  rejection  pulse  is  applied  to  the 
anticoincidence  circuit.  While  the  rejection  pulse  is  present  the  anticoincidence  circuit 
will  not  pass  the  pulses  from  Area  II  and  thus  no  count  can  be  registered  in  the  scaler. 
Therefore,  only  pulses  which  terminate  in  Area  n  are  recorded.  Since  most  Na^'*  and 
Cl^*  traces  terminate  in  Area  II  and  only  a  small  fraction  of  the  background  traces  end 
in  this  area,  the  system  is  selective  with  respect  to  energy. 

The  energy  resolution  for  double  photopeak  traces  on  the  oscilloscope  face  is  approx¬ 
imately  1  percent.  However,  the  energy  gate  at  the  face  of  the  oscilloscope  is  widened 
(Area  II  enlarged)  in  order  to  accept  traces  produced  by  photopeak- Compton  and  Compton- 
Compton  interactions  in  addition  to  the  double  photopeak  interactions.  This  reduces  the 
energy  resolution  but,  at  the  same  time,  increases  the  sensitivity  of  the  system.  In  the 
present  system  the  acceptance  area  (Area  II)  is  approximately  20  percent  of  the  total 
area  within  which  pulses  may  terminate. 

At  the  present  stage  of  equipment  development  it  appears  desirable,  based  on  calcula¬ 
tions  of  the  expected  counting  rate,  to  lower  the  sample  and  the  entire  electronic  apparatus 
in  a  pressure  vessel  to  a  depth  of  several  hundred  feet  in  water  to  greatly  reduce  background 
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countii^  rate  from  cosmic  rays  at  depth  compared  to  the  surface.  This  reduction  in 
background  is  due  to  the  shielding  effect  of  the  sea  water.  The  background  counting  rate 
of  the  system  at  depth  is  0.04  count  per  minute. 

A  multiconductor  cable  is  used  to  supply  power  to  the  equipment  in  the  pressure 
vessel  and  to  allow  a  continuous  readout  of  the  data  on  the  surface.  The  scaler  is  the 
only  piece  of  equipment  shown  in  Fig.  1  which  remains  on  the  surface. 

The  system  is  calibrated  with  the  coincidence  gamma  rays  of  Na^*  (2.75  and  1.38 
Mev)  and  Co®°  (1.33  and  1.17  Mev). 

The  efficiency  of  the  system  for  Na^*  is  0.26  percent.  This  efficiency  is  defined  as 
the  Na^^  counting  rate  in  Area  II  divided  by  the  source  strength.  The  sensitivity  of  this 
system  is  6.9  x  10'^^  curies,  which  corresponds  to  a  signal-to -noise  ratio  of  one. 


EXPERIMENTAL  DATA 

In  June  1960,  this  detection  scheme  was  tested  (2)  with  a  sea-water  sample  from 
beneath  the  Nautilus  while  it  was  docked  at  Portsmouth,  New  Hampshire.  A  tidal  current 
of  about  2  knots  was  flowing  by  the  vessel’s  hull.  A  hose  was  attached  at  the  keel  directly 
under  the  reactor  compartment  in  order  to  collect  the  sea-water  sample.  The  reactor 
was  being  operated  at  50  percent  of  maximum  power.  The  vacuum-tube  prototype  system 
then  in  use  was  not  portable  so  the  sample  was  transported  to  the  U.S.  Naval  Research 
Laboratory  where  it  was  counted.  About  ten  hours  elapsed  between  irradiation  and 
counting.  Essentially  all  of  the  chlorine  had  decayed  by  that  time,  but  the  “signature” 
of  the  radioactive  sodium  was  clearly  visible,  about  50  times  background. 

A  new  transistorized  system  was  developed  in  1961  and  field  tested  at  Key  West, 
Florida,  in  May  1962.  During  these  tests  a  sample  was  taken  from  beneath  the  Nautilus 
while  it  was  docked  at  Key  West,  Florida.  An  estimated  one-quarter  knot  current  was 
flowing  beneath  the  Nautilus.  The  reactor  was  operating  at  10  percent  of  full  power. 

The  5-gallon  sample  was  concentrated  and  then  counted  one  hour  after  extraction  from 
beneath  the  Nautilus.  It  was  not  necessary  to  lower  the  pressure  vessel  as  the  initial 
signal-to -noise  ratio  on  the  surface  was  45. 

The  signal-to -noise  ratio  at  depth  would  have  been  approximately  1000.  A  simple 
calculation  shows  that  the  combined  Na^'*  and  CP*  signal  would  have  been  detectable  on 
the  surface  without  concentration  of  the  salt  with  a  signal-to-noise  ratio  greater  than  two. 

In  a  plot  of  counting  rate  versus  time  the  separate  contributions  of  the  Cl^  and  Na^^ 
isotopes  were  easily  observed.  Detectable  pulses  from  the  Cl^®  isotope  persisted  for 
three  hours  after  extraction  of  the  water  sample  from  beneath  the  Nautilus.  The  Na^*^ 
isotope  was  still  detectable  after  29  hours  with  a  signal-to-noise  ratio  of  6  on  the  surface 
and  150  at  depth. 


DISCUSSION 

Numerous  calculations  (1,3)  have  been  made  of  the  radioactivity  produced  in  sea  water 
by  passage  of  a  nuclear -powered  submarine.  Invariably  the  results  of  these  calculations 
depend  in  a  sensitive  way  upon  many  assumptions,  such  as:  the  total  neutron  leakage  in 
neutrons  per  second  from  the  submarine’s  hull,  the  distribution  of  radioactivity  throughout 


SECRET 


58 


NAVAL  RESEARCH  LABORATORY 


SECRET 


the  wake,  the  volume  of  the  wake  as  a  function  of  time  after  passage  of  the  submarine,  etc- 
The  results  of  most  of  these  calculations  indicate  that  the  present  system  is  probably 
capable  of  detecting  the  nuclear  wake  shortly  after  passage  of  the  submarine. 

Butler  (1)  calculates  the  counting  rate  produced  in  a  system  similar  to  the  present 
one  by  a  sample  removed  from  the  wake  of  a  nuclear  submarine.  It  is  now  possible  to 
improve  the  assumptions  made  by  Butler  and,  in  fact,  to  make  the  calculation  apply  to 
the  present  system.  Using  the  value  recently  determined  by  Miranda  et  al.  (4)  for  the 
total  neutron  leakage  from  the  Nautilus,  and  making  a  set  of  assumptions  which  approxi¬ 
mate  the  conditions  of  the  Key  West  measurement,  one  calculates  a  theoretical  signal-to- 
noise  ratio  of  five  for  the  Key  West  measurement.  The  actual  experimental  value  was 
approximately  1000,  which  is  200  times  larger  than  the  calculated  value.  This  suggests 
that  the  assumptions  made  in  arriving  at  the  calculated  value  may  be  too  conservative, 
and  that  the  signal  being  sought  may  be  more  intense  than  calculations  would  indicate. 

Problems  to  be  solved  include  the  difficulty  in  successfully  extracting  a  sample  from 
the  nuclear  wake.  The  precise  size,  shape,  and  depth  of  the  nuclear  wake  immediately 
aft  of  the  submarine  is  unknown.  At  the  present  stage  of  equipment  development  it  is 
necessary  to  concentrate  the  salt  in  the  sample  by  evaporation  and  to  lower  the  concen¬ 
trated  sample  in  a  pressure  vessel  to  reduce  the  background  counting  rate  by  a  factor  of 
approximately  23.  Since  these  two  operations  require  about  one -half  hour,  one  must  probe 
the  wake  for  a  sample  and  later  process  it  without  any  assurance  that  the  sample  was 
actually  extracted  from  the  nuclear  wake.  The  delayed  readout  of  this  prototype  model 
of  the  detector  system  greatly  increases  the  difficulty  in  obtaining  a  sample  from  the 
nuclear  wake.  It  is  hoped  that  the  sensitivity  of  future  prototype  systems  will  be  suffi¬ 
ciently  increased  so  that  it  will  not  be  necessary  to  concentrate  the  sample  or  lower  it 
in  a  pressure  vessel. 

The  practical  application  of  this  system  to  the  detection  of  a  nuclear  wake  depends 
on  the  spatial  distribution  of  the  wake  as  a  function  of  time.  There  is  some  experimental 
evidence  that  the  nuclear  wake  remains  a  well-defined  entity  for  several  hours  after 
passage  of  a  submarine  (5). 

The  ability  of  this  system  to  detect  a  nuclear -powered  submarine  is  directly  dependent 
on  the  neutron  leakage  from  this  submarine.  In  the  design  of  this  system  it  has  been 
assumed  that  the  neutron  leakage  of  a  potential  enemy  submarine  is  the  same  as  from 
our  submarines. 

A  simple  modification  in  the  present  system  could  enable  it  to  count  the  Cl^*  and  Na^'* 
gamma  rays  separately.  From  the  relative  counting  rates  of  the  two  isotopes  one  can 
calculate  the  elapsed  time  since  passage  of  the  nuclear  submarine.  This  modification  has 
not  been  made  in  the  present  system  because  the  main  problem,  at  present,  is  detection 
of  the  submarine. 

Although  many  problems  exist,  particularly  in  successfully  extracting  a  sample  from 
the  nuclear  wake,  these  problems  do  not,  at  present,  appear  insurmountable  and  work  is 
continuing.  The  next  step  will  involve  extracting  a  sea- water  sample  from  beneath  the 
reactor  of  a  nuclear  submarine  while  it  is  underway.  When  this  step  is  successfully 
concluded,  wake-sampling  equipment  will  be  developed  to  be  operated  from  a  surface 
ship,  and  samples  will  be  taken  from  the  nuclear  wake  aft  of  the  submarine. 

Many  improvements  can  and  will  be  made  in  this  equipment.  It  is  hoped  that  it  will 
prove  to  be  valuable  not  only  in  classifying  nuclear  submarines,  but  also  in  providing  a 
new  technique  for  the  study  of  the  diffusion,  drifting,  and  dissipation  of  wakes. 
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